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Effect of Packing Configuration 


on Mass and Heat Transfer 


in Beds of Stacked Spheres 


L. R. GALLOWAY?, 


Experimental measurements are reported on the simultaneous rates 
of mass and heat transfer for the surface evaporation of water into air 
from two orthorhombic, a rhombohedral, and two simple cubic assem- 
blages of uniform celite spheres, in the particle Reynolds number range 
300-1200. The measurements were confined to the pre-determined 
constant rate drying period of the spheres. 

No appreciable differences in mass or heat transfer factors were 
found for the various assemblages, except for the cubics, which exhibited 
slightly lower transfer factors than the other assemblages at the higher 
Reynolds numbers. The negative slope of the log j vs. log Re curves was 
designated the “index of turbulence inhibition”. This index was cor- 
related with the projected fractional free area of the packings, as well 
as with the wall effect. 

Various methods proposed by investigators of randomly packed beds 
for modifying the mass transfer factor and/or the Rzynolds number 
failed to improve the overall correlation of the results. The inapplic- 
ability of a simple analogy between momentum and mass transfer, as 
proposed by Ergun, is demonstrated. 


Based on the model of McHenry and Wilhelm, the magnitude of 
the axial mixing factor required to correct for the use of logarithmic 
mean potentials is indicated. 


oo study is concerned with the 
simultaneous transfer of mass and 
heat between a fluid flowing through 
a packed bed of spheres and the sur- 
face of these spheres. It does not con- 
cern transfer to or from the walls 
enclosing the bed. 


Important as they are in such oper- 
ations as drying, leaching, pebble 
heating, adsorption and ion exchange, 
as well as in heterogeneous reactor 
design, packed beds of perfect or 
nearly perfect spheres undergoing an 
exchange of mass and/or heat with 
fluid permeating the void spaces be- 
tween the spheres have been experi- 


mentally investigated before (-1*), 
The studies which have been re- 
ported, however, all deal with ran- 
domly packed and hence inherently 
unreproducible packing assemblages. 
In the present investigation a new 
variable has been introduced, namely, 
the geometric configuration of the 
packing. This was done by dealing 
with carefully stacked, rather than 
randomly packed, assemblages. 


Fluid Friction 


A precedent for this type of study 
was set by Martin et al. (9) who 
investigated the fluid friction charac- 
teristics of nine different orderly 
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configurations of brass balls, stacked 
according to the arrangements elabo- 
rated by Graton and Fraser ‘2°), Four 
of these arrangements, illustrated 
schematically in Figure 1, were 
studied in the present work. The 
simple cubic represents the loosest, 
while the rhombohedral (hexagonal 
close-packed) represents one of the 
two densest, stable arrangements of 
uniformly sized spheres. The two 
orthorhombics have identical lattices, 
but differ by 90° in their orientation 
with respect to the major flow axis. 
Thus the plan of No. 2 Clear Passage 
is equivalent to the elevation of No. 
4, while the elevation of No. 2 Clear 
Passage is equivalent to the plan of 
No. 4. Orthorhombic No. 2 Blocked 
Passage was not constructed in the 
present work. 


Martin’s results on the four assem- 
blages, in the particle Reynolds num- 
ber range of the present investigation, 
are shown in Figure 2. While the 
large differences in friction factor 
between the cubic, orthorhombic and 
rhombohedral packings can be ex- 
plained in part at least by the large 
fractional void volume | differences 
between these three arrangements, the 
four-fold difference in friction factor 
between the two orthorhombics is a 
surprising result in view of the 
identical lattice of the two packings. 
This result was confirmed in the 
present work by making longitudinal 
pressure profiles on the two ortho- 
rhombic assemblages in conjunction 
with a number of the mass and heat 
transfer measurements subsequently 
reported. The data obtained for air 
flow through celite spheres are plot- 
ted in Figure 2, where they show 
good agreement with Martin’s best 
lines for liquid flow through metal 
spheres. 
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Figure 1 — Arrangements of spheres 
studied. 
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Whether the inordinate increase in 
friction factor from Orthorhombic 
No. 2 to Orthorhombic No. 4, which 
represents a corresponding increase 
in rate of momentum transfer, would 
reflect itself, in part at least, in an 
increase in rates of heat and mass 
transfer between the spheres and the 
fluid, was one of the questions which 
spurred the present study. It should 
be noted that this large percentage 
increase in friction factor cannot 
simply be attributed to form drag, 
since it was shown by Martin to be 
even larger in the viscous flow range, 
where skin friction alone prevails. 
This is in contradistinction to the case 
of No. 2 Clear Passage versus No. 2 
Blocked Passage, where a consider- 
ably smaller difference in friction 
factor occurs only at particle Rey- 
nolds numbers exceeding 20 and the 
difference increases with Reynolds 
number. 


Experimental Method 


The spheres used to construct the 
packings were composed of 50% 
celite, 37.5% alundum cement and 


Martin—Rhombohedral No.3 


y—-Martin— Cubic No.| 


@ this study -Orthorhombic No. 2 


@ this study—Orthorhombic No. 4 
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Figure 2—Fluid friction characteristics of packings. 
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12.5% colloidal kaolin, by weight. 
The method developed for fabrica- 
ting them has already been reported 
(21), They had a mean diameter of 
0.673 in., with a standard deviation of 
less than 1° from this mean. Their 
high internal porosity allowed them 
to retain about two-thirds their own 
volume of water, which was trans- 
ferred to an air stream during the 
course of an experimental run. The 
spheres also fulfilled the requirements 
of being tough, machineable, and 


subject to negligible attrition — in 
handling. 
The assemblages were built by 


individually drilling small holes in 
each sphere and pinning each to its 
tangent neighbour with 0.022 in. 
stainless steel fishing wire. The wire 
was secured in each hole with Aral- 
dite cement. The tendency of the 
spheres to bridge and cause extra 
porosity at the wall was eliminated by 
using fractional spheres at the walls 
and corners, in the manner of Martin 
et al .{'%, The bundles of spheres 
were enclosed on all sides except the 
top and bottom by 1/16-in. brass 
plate cemented to the faces of the 
fractional spheres. Photographs of 
the assembled packings are shown in 
Figure 3. 

Packing Nos. 1, 2, and 3, based as 
they are on square horizontal layers 
(see Figure 1), were housed in square 
containers with half-spheres at ‘the 
walls and quarter-spheres in the cor- 
ners, while Orthorhombic No. 4, 
composed of rhombic layers, could 
only be adapted to a hexagonal con- 
tainer, using half-spheres at the walls 
and one-third spheres in the corners. 
A fifth stacked assemblage, Cubic No. 
1A, was also constructed and studied. 
This was similar in all respects to 
Cubic No. | except that no fractional 
spheres were used at the walls. The 
effect of exposing the additional wall 
area to the air flow could thus be 
ascertained. 

Dimensions of the various assem- 
blages are shown in Table 1. Eight 
layers of spheres were used for each 
stacked assemblage, a compromise 
between too few layers which could 
make for end effects and too many 
layers which would cause thermal 
and diffusion equilibrium to occur 
part way through the packing, thus 
rendering the determination of heat 
and mass transfer coefficients by 
integral methods impossible. The var- 
lation in fractional void volume was 
than double that which can 
normally be obtained in randomly 
packed beds of uniform-size spheres 
(22), Measurements were made 


more 
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| on the spheres randomly packed, to 

ensure that the methods used give 

re resu'ts in reasonable agreement with 

of other investigators of transfer in 
of randomly packed beds. 

The apparatus in which the pack- 
ings were housed is illustrated sche- 
matically in Figure 4. After removing 

any entrained moisture from the 
pressure- regulated building compres- 
re | sor air, this air was throttled to 
reduce its relative humidity to below 
40°, orifice-metered into the test 
| column, and finally exhausted to the 
atmosphere. Both the square and the 
hexagonal test columns were made 
| from 4-in. aluminum plate, with 4-in. 
aluminum flanges. The inside cross- 
sectional dimensions of these columns 
were such as to make a snug fit when 
the assemblages of spheres with brass 
1e side plates were placed inside the 
ra columns. Empty calming zones on 
either side of the assemblages, and 
about twice their depth, were lined 
with 1/16-in. brass plates inorder to 
maintain a constant cross-section 
throughout the entire column. Five 
equally spaced pressure taps (see 
Figure 3f) from each assemblage 
were brought through a flange of the 
column via a brass plate attached to 
the wall containing the pressure taps. 
These were connected to open and 
differential draft gauge manometers, 
in order to measure in conjunction 
with barometer readings, the absolute 
pressure and its variation throughout 
the depth of packing. The columns 
were insulated with two inches of 

glass wool. 


The dew point of the entrance air 
and of the air after it had picked up 
moisture from the wet spheres, was 
measured to within 0.3 Fahrenheit 
degrees by passing samples of these 
streams through a Foxboro “dewcel” 
linked to a continuous recorder. 
Temperatures immediately below and 
immediately above the packing were 
measured with carefully calibrated 
thermometers, reading to the nearest 
0.1 Fahrenheit degree. 


Before making a run the entire 
packing assemblage was soaked in a 
water bath held within 2 Fahrenheit 
degrees of the adiabatic saturation 
temperature of the inlet air, for at 
least three hours. After shaking off 
excess water, the wet assemblage was 
inserted into the column and the run 
started. All the relevant measure- 
ments were made during the constant Figure 3—Photographs of assemblages 
rate period of drying. This was pre- 

5 determined at several air velocities (a) eubie no. | (d) rhombohedral no. 3 
l through two layers of randomly (b) cubie no. 1A (e) orthorhombic no. 4 
packed spheres, a typical drying (ec) orthorhombic no. 2 clear passage (f) pressure taps 

; curve of which is illustrated in Figure 
5. Even at the highest particle Rey- 
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TABLE 1 
CHARACTERISTICS OF 


PACKINGS 







Side or Cross- Surface 
Packing Shape of Diameter | Sectional No. of Area of | Depth 
Arrangement Column of Column Area Spheres Spheres inches 
inches sq. ft. sq. ft. 
Cubic No. 1 & 1A Square 4 11/16 0.1526 392 3.873 5.38 
Orthorhombic No. 2 2 7 
clear passage Square 4 11/16 0.1526 392 3.873 +.66 
Rhombohedral No. 3 Square 4 11/16 0.1526 392 3.873 3.81 
Orthorhombic No. 4 Hexagonal 2 11/16 0.1303 384 3.795 5.38 
Random Circular 3.00 0.0491 142 1.404 
Random Square 4 11/16 0.1526 142 | 1.404 
Draft gauge measuring jee Ae 
pressure drop through packing. a ‘ na 


4p 20 

[ | Orifice | 
Upstream pressure | | | 
| 
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Orifice 


“ 
——— =— From ‘Dewcel. 


k- ies 


| 


Pressure at bottom 
of packing 
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~* To atmosphere 





Thermometer 
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Thermometer 


Entrainment ly. t 
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Pressure 
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Figure 4—Schematic of apparatus. 


nolds number studied, about 1,200, 
this period lasted more than one-half 
hour, giving ample time for reliable 
measurements. 

Most of the runs performed were 
with unheated air, but for a few runs 
the inlet air had to be pre-heated 
order to reduce its relative humidity. 
The outlet air was thus kept at a 
relative humidity below 949%, usually 
in the range of 80-90% 


Data Processing 

The following rate equations were 
calculate the 
mass and heat transfer respectively: 


used to coefficients of 


N = kgA(p*—p)1.m (1) 
q hA(t —t*)im (2) 
Phe quantity p* was taken as the 


at the tem- 
‘The mean partial pressure 
derived by Gamson et al. 
simplifies to the logarithmic 
mean used in equation | when p*, Py 


vapor pressure of water 
perature t* 
difference 
(4) 
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and p, are small compared to the total 
pressure P as in the present work, 
where the two means never differed 
by more than 0.4%. The partial pres- 
sures of the water vapor in the inlet 
and outlet streams were taken as the 
saturated vapor pressures correspond- 
ing to the respective dew points. In 
conjunction with the measured abso- 
lute pressures of these streams, the 
partial pressures were used to estab- 
lish the humidities. The mass transfer 
rate N was evaluated from. the 
measured air rate and the change in 
humidity, while the heat transfer 
rate q was evaluated from the cor- 
responding thermal data. 

The transfer area A was taken as 
the geometric surface area of the 
spheres. No correction was applied 
for the area lost due to insertion of 
the connection wires, since this loss 
amounted to only 1% 


The quantity t*, the surface tem- 


perature of the spheres during the 
constant rate dryi ing period, was as- 
sumed constant throughout the depth 


Largest Smallest Projected 

Fractional Fractional Fractional Fractional 
Void Free Free Free 
Volume \rea \rea \rea 
0.4764 1.0 0.215 0.215 
0.3954 0.635 0.215 0.043 

0.2595 0.349 0.215 0 

0.3954 1.0 0.093 0.093 


of packing and equal to the adiabatic 
saturation temperature of the outlet 
air. This calculated from the 
measured temperature, dew — point, 
and absolute pressure of the exit air 
and fell within 0.5 Fahrenheit degrees 
of the adiabatic saturation 
ature calculated from the 
data. It also fell within 0.5 Fahren- 
heit degrees of the values measured, 
for several trials, by baking thermo- 
couples into the spheres with junc- 
tions located at the surfaces opposite 
to the entrance points of the insulated 
thermocouple wires. The same as- 
sumption was made by Gamson, 
Thodos and Hougen ‘*) in their sim- 
ilar measurements on randomly packed 
celite spheres and cylinders in a 
particle Reynolds number range en- 
compassing and exceeding that ‘of the 
present study, as well as by Wilke 
and Hougen ‘?3) on celite cylinders 


was 


temper- 
inlet air 


at lower Reynolds numbers. The 
latter investigators resorted to this 
assumption because of difficulties in 
making accurate measurements of 
surface temperature by means of 
attached thermocouples. In the dis- 
cussion of the Gamson-Thodos- 
Hougen paper it was noted by 
Chilton that Hurt’s experimental 


check of t* “with the temperature of 
adiabatic saturation, or the wet bulb 


temperature, was as good as the 
agreement is between these two 


temperatures” Although Hobson and 
Thodos ‘') subsequently criticized 
this assumption for very low Reynolds 
numbers, they found that it was 
increasingly applicable for the air- 
water system as the Reynolds number 
approached the lower values of the 
present investigation, even under 
non-adiabatic drying conditions. 
More recently Heertjes and Ringens 
(24) have shown experimentally that 
the constant rate evaporation of 
four different liquids from earthen- 
ware blocks into an air stream, cover- 


for 
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ing wide range of air velocity, the 
solid temperature was equal to the 
wet ulb temperature. It should be 
note! that in the present study the 
wet-bulb temperature is interchange- 
able with the adiabatic saturation 
temperature. This follows from the 
fact that thermal radiation to the 
pack ing surface was almost wholly 
absent, since most of the spheres 
could “see” nothing but other spheres 
at the same surface temperature. If 
the value of the air-water psychro- 
metric ratio is taken as 0.236, recently 
reported as the most reliable value of 
this ratio in the absence of radiation 
effects (*5), then for the low humidity 
range of the present study the wet- 
bulb temperature falls within 0.1 
Fahrenheit degree of the adiabatic 
saturation temperature. 

The data for each packing were 
correlated by means of the particle 
Reynolds number and the j-factors of 
Colburn (6) and Chilton (27), Phys- 
ical properties of air-water vapor 
mixtures were obtained from the most 
reliable sources available (4, 28, 29, 30), 
For all the runs performed the 
Schmidt and Prandtl numbers were 
within 1% of their average values of 
0.609 and 0.718 respectively. 


Experimental Results and Their 
Interpretation 

The results for 57 runs on Assem- 
blage Nos. 1, 2, 3, and 4 are summar- 
ized in Figures 6 and 7. Except for 
the lower than average points of 
Cubic No. 1 at the higher Reynolds 
numbers, the points for the four 
assemblages largely overlap. The 
straight lines drawn in Figures 6 and 
7 have been obtained by fitting the 
data for each assemblage, using the 
method of least squares, by equations 
of the form 


ja = aRe™ (3) 
jn = bRe™ (4) 


and 


The empirical equations for each as- 
semblage along with the correspond- 
ing standard errors are recorded in 
Table 2. Even for Assemblage Nos. 
1-4 taken together the standard errors 
in jg and j, were + 15% and + 18% 
respectively. It is therefore concluded 
that packing configuration plays a 
much smaller role in mass and heat 
transfer than it does in fluid friction. 

The negative slopes of the log j vs. 
log Re curves, m and n in equations 
3 and 4 respectively, are an inverse 
measure of the development of 
transfer-effective turbulence in a flow 
system, varying as they do from 
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Figure 5—Drying curve for two layers of randomly packed spheres at Re = 443. 
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Figure 6—Effect of packing arrangement on mass transfer. 


unity for viscous flow to zero for 
the rare case in packed beds of fully 
developed turbulence. Coppage and 
London (6) were able to correlate n 
with ¢ for screen matrices, broken 
solids, and random beds of spheres. In 
the present experiments, however, of 
the interstitial properties listed in the 
last four columns of Table 1, the 
projected fractional free area gave 
the best correlation with these slopes. 
This correlation is shown in Figure 8, 
where the arithmetic mean of m and 
n is plotted against the projected 
fractional free area for Assemblage 
Nos. 1 and 4, and the corresponding 
point for an isolated sphere is also 
included. The average slope for an 
isolated sphere was calculated for 
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Re = 300—1200, Sc = 0.609, from the 
Fréssling-Ranz equation @), The 
ordinate of Figure 8 may be con- 
sidered an “index of turbulence inhi- 
bition”. Its steady increase with the 
area of unobstructed passageways in 
the net direction of fluid flow sug- 
gests that channeling occurs through 
these passageways at the expense of 
cross-flow. The rhombohedral _ar- 
rangement with its completely 
blocked off “view” represents one 
extreme with a maximum of tortu- 
osity and therefore a maximum de- 
velopment of turbulence, while the 
isolated sphere with its surrounding 
infinite flow area represents the other 
extreme with a minimum develop- 
ment of turbulence, for the given 
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TABLE 2 
EMPIRICAL EQUATIONS FOR PACKINGS 
Packing Mass Standard Heat | Standard 
Arrangement Transfer Error Transfer Error 
Cubic No. 1 ja = 0.460 Re-®-*49 | + 3.6% | jn = 0.612 Re-®:379 | + 7.9% 
Orthorhombic No. 2 ja = 0.084 Re~*-925 | + 9.0% | jn = 0.097 Re-°-980 | 412.7% 
clear passage 
Rhombohedral No. 3 ja = 0.064 Re*9-9°S | + 8.6% | jn = 0.088 Re~°-°23 | 414.0% 
Orthorhombic No. 4 ja = 0.130 Re-®-""4* | + 6.5% | jn = 0.279 Re-9-213 | + 7.9% 
Combined No. 1,2,3,4 ja = 0.135 Re~®-"7 | 414.6% | jn = 0.197 Re--189 | 418.3% 
Cubic No. 1A ja = 2.963 Re~°-638 | + 9.2% | jn = 5.423 Re-9-722| 4+ 8.9% 


particle Reynolds number range and 
in the absence of wall effect. 


” 


The “turbulence inhibition index 
of deep randomly packed beds of 
spheres in the particle Reynolds 
number range of the present study is 
0.30 or very close to this value, 
according to several investigations 
5. 6, 7, 9, 11, 12, 16, $2), This value is 
bracketed by those for Orthorhombic 
No. 4 and Cubic No. 1. Since the 
void fraction of random beds of uni- 
formly sized spheres is usually be- 
tween that for the orthorhombic ar- 
rangement and that for the cubic (?), 
it can be surmised that a heat and 
mass transfer model for such beds is 
a mixture of local cubics and local 
orthorhombics. 


Because slightly better correlation 
was obtained for }, than for },, subse- 
quent discussion will be largely re 
stricted to mass transfer, though it 
should be realized that similar state- 


ments can be made about heat trans 
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fer. Sherwood (33), in fact, has 
pointed out that taking the particle 
surface temperature as the adiabatic 
saturation temper rature is in effect 
saying that the ratio of }, to jq is 
given by (Pr/Sc)?/%. This is borne 
out by the fact that while (Pr/Sc)?/3 
for the present study is (0.718/0.609 )?/3 
or 1,120, the average value of },/jq 


for the 57 runs discussed was 1.128. 


Figure 9 compares the best straight 
line for the combined results of 
Assemblage Nos. 1-4 with corre- 
sponding results on randomly packed 
beds. As representative of the latter, 
the line based on the empirical 
equation of Gamson, Thodos and 
particle Reynolds 
numbers exceeding 350 was chosen, 
since investigators also dealt 
with celite spheres and the air-water 


Hougen ‘for 
these 


system. The two lines cross each 
other and deviate on the average by 
only about 10% over the range in- 


volved. Moreover, the relatively low 
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Figure 8—Dependence of turbulence 
inhibition on passageways open to 
unobstructed through-flow. 


negative slope of the line for the 57 
runs on the four assemblages can be 
explained by the non-selective aver- 
aging method used. If, for instance, 
more weight had been given to the 
cubic assemblage, the fractional void 
volume of which is very close to that 
of “random-loose” beds (22), and less 
to the rhombohedral, the void frac- 
tion of which is more than 0.11 below 
that of even “random-dense” beds (22), 
then this slope would have markedly 
increased. On the other hand, the 
slope of the Gamson-Thodos-Hougen 
line is somewhat higher than that ob- 
tained by investigators of other 
systems involving randomly packed 
spheres in the same or even lower 
particle Reynolds number range 
is 6. 3s 0 B. . 36s 16), This has been 
explained (71) by the shallowness of 
the packed beds used by Gamson et 
al.—as low as three layers of spheres 
in some cases—an effect which pre- 
vents the maximal development of 
turbulence. 


The dozen runs on randomly 
packed beds made in the present 
study are also summarized in Figure 

The upper six points, obtained 1 
the square column, are in good agree- 
ment with the best line of Gamson 


et al. The six points obtained in the 
3-in. diameter round column fall 
somewhat lower, though they are 


still above the line for an_ isolated 
sphere by a factor of two. The low 
points are readily explained by the 
small ratio of column diameter to 
sphere diameter for this case, that is, 
by the large wall effect. The marked 
tendency of the spheres to bridge at 
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the wall causes excessive voids and 
hence the channeling of an undue 
roportion of the air near the wall 
(8, 8%, 35), This channeling air cannot 
in large measure participate in the 
exchange of mass and heat with the 
spheres, thus lowering the integral 
coeflicients of mass and heat transfer. 
The same effect has been found for 
heat transfer in the absence of mass 
transfer (% 11), 

Figure 10 compares the results on 
Cubic No. 1A, which makes only 
point contact at the walls, with those 
for Cubic No. 1, which contains the 
flat side of half-spheres in contact 
with the walls and thus removes 7/4 
or 78.5% of the wall area. An even 
greater negative slope ' obtained 
for No. 1A than for No. 1. This can 
be explained by the alae inhib- 
iting action of the added wall area, 
which is in the direction of the net 
air flow and therefore adds to the 
skin friction. The two lines intersect 
at a particle Reynolds number close 
to that at which Rose and Rizk (35) 
found that the wall-effect factor for 
fluid friction switched from positive 
to negative. A similar intersection 
due to wall effect on j,, though at a 
higher value of Re, was found by 
Denton, Robinson, and Tibbs in their 
differential measurements on steady 
state heat transfer from electrically 
heated copper spheres in packed 
beds (7% 11), A positive deviation from 
the case where wall effects are absent 
denotes the predominance of the wall 
surface resistance over channeling 
effects, while a negative deviation is 
attributable to the predominance of 
channeling. Although such channeling 
occurs along the wall in the case of 
a randomly packed bed with a small 
value of D./D,, it probably occurs 
through the centre of the bed in the 
case of Cubic No. 1A, in which the 
void fraction near the wall is no 
larger than in the core of the bed ‘*), 


Correlating Functions 

An attempt was made to narrow the 
spread of the points in Figure 6 by 
means of modifications to Re and Ja 
which have achieved some success in 
the case of randomly packed beds. 
These are listed chronologically in 
Table 3. Most of the modifications 
were derived on the assumption that 
the effective velocity to be used in Re 
and sometimes also in j, is the average 
interstitial velocity V,/e rather than 
the superficial velocity V,. In some 
of these cases six times the hydraulic 
diameter, D,/(1-e), also replaces the 
particle diameter D,, in the Reynolds 
number. Ranz (4!) used the maximum 


0.100 


0.080 


0.060 
® 


~ 0.040 


ee 


0.030 


0.020 


Gamson et al 


Asolated Sphere 


(Sc = 0.609) 


O Random, square column, De/Dp = 7.0 


@ Random, round column, De /Dp #45 


300 400 


Figure 9—S:acked bed composite, 





600 800 
Re 


random beds, and isolated sphere mass 


1000 1200 1500 


transfer compared. 





0-10 
0-09 


0-08 


0:04 








—O— Cubic No.1 ‘@ | 
» *< | 
_— -@ = | 
0:03 Cubic No.ia | 
| 
0:025 - 1 - ee 1” aa ! —_ | = al 4 
100 150 200 300 400 $00 600 800 1000 1500 


Figure 10—Effect of wall on mass transfer. 


interstitial velocity V,/a in his modi- 
fied moduli. The use of the Schmidt 
number raised to some power to 
modify j, merely signifies that the 
given investigators employ a power 
for Se other than the 4 originally pro- 
posed by Chilton and Colburn. Its 
additional use in the modified Rey- 
nolds number of Ranz arises out of 
the form of the Frdéssling equation 
for mass transfer from an_ isolated 
sphere (7), These Schmidt number 
factors may be ignored in attempting 
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to correlate the present data, which 
were all at constant Sc. 

Figure 11 shows the data plotted 
according to the method used by 
McCune and Wilhelm to bring their 
liquid-solid mass transfer data on 
fixed beds into line with those for a 
single sphere (6). No success is 
either in reducing the 
scatter of the points of Figure 6 or in 
aligning the points about the isolated 
sphere line. 

Phe method equivalent to that used 


achieved 
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TABLE 
DIMENSIONLESS MODULI FOR Mass TRANSFER IN RANDOMLY PAacKEL p BEDs OF SPHERES 


Investigators 
McCune : ond Wilhelm (6 
Gaffney and Drew (9 
Gamson (38) 
Ergun (32 
Ranz (31) 


Chu et al. (13, 45 


© Cubic No.! i 

@ Orthorhombic No.2 is 
— @ Rhombohedral No.3 . 
@ Orthorhombic No.4 ® 


[S) 


2.015 -—- 
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Year Modified Re Modified ja 
1949 Re € i — 
1950 Re/e ja Sc~°-0 
1951 Re/1—e ja/(1 —€)° 
1952 Re/1—e jae Sc*/8 
1952 Sc?/3 Re/a jaa Sco!/8 
1953 Re/1—e ja 

—— 
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Figure 11—Correlation after McCune and Wilhelm. 
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Figure 12—Correlation after Gaffney and Drew. 


by Gaffney and Drew ‘®) to correlate 
liquid- -solid mass transfer data in 
randomly packed beds of spheres is 
employed in Figure 12. No rectifica- 
tion of the data results. The straight 
line drawn is based on an empirical 
equation derived by Ergun ‘%) to 
represent the plot made by Gaffney 
and Drew of their own data, as well 
as those of McCune and Wilhelm. It 
misses the present data completely. 
Figure 13 illustrates the application 
of Gamson’s (88) modified moduli to 
the stacked bed results. The arbitrarily 
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empirical use of the factor (1-¢)®? in 
the denominator of the ordinate, 
especially after I-¢ is already in the 
denominator of the abscissa, detracts 
from the rationality and _ sensitivity 
of this method. Little or no rectifica- 
tion of the data is achieved over that 
already present in Figure 6. More- 
over, the line based on Gamson’s 
correlation equation for fixed and 
fluidized beds of spheres is inade- 
quately representative of the present 
data. 

The combined use of interstitial 


velocity and hydraulic diameter in 
the dimensionless groups, as proposed 
by Ergun (82), is the mass transfer 
analogue of the method developed by 
Blake (9), Kozeny ‘#), Fair and 
Hatch (41), Carman ‘#) and Ergun 
(48) for fluid friction in granular 
media. That it impairs rather than 
improves the correlation of the pres- 
ent data is apparent on comparing 
Figure 14 with Figure 6. In addition, 
the mass transfer equation derived by 
Ergun, equation 5b below, fails by 
an average factor of about 16 to pre- 
dict the present data. Even if perfect 
mixing in the beds is assumed, an 
extreme possibility suggested by 
Ergun for gas flow through shallow 
beds, his line still lies at an aver age 
factor of approximately 5 above the 
recalculated points, which show 
more scatter than ever. 


The modified dimensionless moduli 
proposed by Ranz (31), like those of 
McCune and W ithelm, are designed 
to reduce to their corresponding un- 
modified forms for the case of an 
isolated sphere, for which « = a = 1, 
The present data are plotted accord- 
ing to this scheme in Figure 15, 
which again shows no improvement 
over the spread in Figure 6, but does 
have the virtue of approximately 
aligning the points about the isolated 
sphere line. It would seem then that 
for the stacked beds of this study, a 
is a more successful parameter than ¢ 
for relating the mass transfer data to 
those for an isolated sphere. The 
relevant difference between stacked 
and randomly packed beds for trans- 
fer processes is probably the contin- 
uous but periodic variation of frac- 
tional free area with bed height in 
the case of the former, as opposed to 
the statistical constancy of fractional 
free area in the case of the latter (44). 
A constant fractional free area is 
moreover equal in value to both a 
and ¢. This probably accounts for the 
success of McCune and Wilhelm in 
correlating their randomly packed 
bed data with that for a single sphere 
using ¢« as parameter, despite the 
failure of their method in the present 
instance. 

Finally, the equation of Chu et al. 
(18. 49), which has provided a general- 
ized correlation of mass transfer in 
fixed and fluidized beds of granular 
particles, is tested in Figure 16. Here 
again there is no reduction in the 
scatter of the points, and Chu’s line, 
like Gamson’s, falls somewhat high. 


It is concluded, therefore, that 
methods appropriate to the correla- 
tion of mass transfer data on ran- 
domly packed beds are inappropriate 
to stacked beds. The same conclusion 
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for fluid friction was arrived at by 
Maitin (194), 


Mo.nentum-Mass Transfer Analogy 
Starting with Reynolds’ original 
analogy (45), Ergun (82) has derived a 
“theoretical” relationship between 
momentum and mass transfer in 
packed beds, ignoring completely the 
fact that the Reynolds and derivative 
analogies apply only to skin friction, 
but not to form drag (47), Ergun’s 
analogy may be expressed as 
«3 
6eSc"-tjq = 26 (5) 
1—e 
Liquid-solid data ‘*-*. 9) were used to 
evaluate the index n as unity. The 
simplified relationship is then 


ee eee (Sa) 
3 1—e 
Combining equation 5a with his well- 
supported equation for fluid friction 
in randomly packed beds (4%), Ergun 
finally arrived at 

Sk 150 (1—e) ss - 

6€Sc!/3jq = + 1.75 (5b) 
Re 

When equation 5b was applied to 
gas-solid data (4 23, 48,49), it failed 
completely to predict the experi- 
mental results by factors ranging 
from 5 to 30. The smallest factor in 
this range applied only when the 
extreme assumption of perfect mix- 
ing was employed for the shallow 

bed mass transfer data of Hurt (48), 


On a purely empirical basis, Chu et 
al. (1% 45) have proposed “for the 
practical purpose of engineering de- 
sign” the following approximate 
analogy between mass and momentum 
transfer in fixed and fluidized beds of 
granular solids: 


mane (6) 
"10 1-€ 
Both equations 5a and 6 imply that 
ja is directly proportional to f if the 
other terms in the equations are fixed. 
This is completely contradicted by 
the results on the two orthorhombic 
assemblages. These have the same 
void fractions, and since Sc was 
essentially constant for the present 
experiments, all the terms other than 
j, and f were thus fixed. Nevertheless 
a four-fold increase in f from Ortho- 
rhombic No. 2 to Orthorhombic No. 
4 was accompanied by no increase in 
Jas 
The application of equations 5a 
and 6 to the present assemblages at a 
particle Reynolds number of 1,000 is 
shown in Table 4. While the E rgun 
analogy yields values of j, which are 
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assumption of perfect mixing). 
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TABLE 4 
TEST OF MOMENTUM-Mass TRANSFER ANALOGIES AT Re = 1000 


f 


Packing Experimental 


Ergun analogy 


ja ja _ 
Chu analogy Experimental 


















Arrangement Martin et al. Equation (5a) | Equation (6) Table 2 
(19) Sc = 0.609 
Cubic No. 1 1.8 0.307 0.0373 0.0439 
Orthorhombic No. 2 
clear passage 5.8 0.590 0.0593 0.0706 
Rhombohedral No. 3 30 1.07 0.0708 0.0677 
Orthorhombic No. 4 23 2.34 0.235 0.0592 
= 7 7 7 
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Figure 15—Correlation after Ranz. 


high by factors ranging from 7 for 
Cubic No. 1 to 40 for Orthorhombic 
No. 4, the Chu analogy predicts to 
within 16% or less the values of j, 
for the first three assemblages, though 
of course it still misses the j, for No. 
4 by a factor of 4. 


The inapplicability of Ergun’s 
analogy to packed beds can in part at 
least be explained by its failure to 
differentiate between skin friction 
and form drag, only the former of 
which has appreciable effect on mass 
and heat transfer (5), It is known 
that for flow past single cylinders (51), 
single spheres (47) and closely array od 
cylinders ‘: 5), turbulent wake for- 
mation results in a much larger dissi- 
pation of energy by 
than by skin friction, 
Reynolds numbers. 
follow then, 


form resistance 

at the higher 
It would seem to 
that for closely arrayed 
or packed beds of spheres at com- 
parable Reynolds numbers, turbulent 
wakes form behind the spheres, caus- 
ing a large dissipation of energy (high 
fluid friction) without a correspond- 
ing increase in mass or heat transfer. 
This would 
that 


also account for the fact 
eventually becomes 


independent of Re _ for 


whereas f 
randomly 
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packed beds (7:!!: 4), | continues to 
fal] (1-18, 23, 48, 49) right up to Re = 
60,000, the highest particle Reynolds 
number at which mass or_ heat 
transfer has been studied ‘7 11), The 
same effect again holds true for flow 
past isolated spheres and cylinders, as 
well as for flow through rough pipes 
(47, 50,53), Flow through smooth 
pipes, by contrast, displays only skin 
friction and therefore both f and j 
continue to fall parallelwise up to the 
highest of measured pipe Reynolds 
numbers (47), as do the calculated 
skin friction factor (5!) and experi- 
mental mass and heat transfer factors 
for flow normal to single cylinders (#7). 
KE xponential equations of the form of 
equations 3 and 4 are therefore em- 
pirically more applicable in the 
transitional and turbulent range of 
particle Reynolds number than the 
type similar to equation 5b. 


Mixing 

The common practice of using 
logarithmic mean driving forces in 
evaluating transfer coefficients, as in 
equations 1 and 2 of the present 
study, implies the absence of any 
longitudinal mixing in packed beds. 
Ergun ‘3*) has suggested that the 


Assumptions 


opposite extreme of perfect mixing 
is a better assumption for gas flow 
through shallow beds, which would 
require basing the transfer coefficients 
on exit driving force. However, his 
analysis of Hurt’s data (48) in the 
transitional range of Reynolds num- 
ber which led him to this conclusion 
was based on the form of equation 
5b, which as indicated above, is less 
applicable than equation A similar 
analysis of Hurt’s data based on 
equation 3 leads to the opposite con- 
clusion, as shown by the straight lines 
obtained by Resnick and White (9) 
on plotting these data as log j, versus 
log Re, basing j, on logarithmic mean 
driving force. 

Figure 17 shows the mass transfer 
data on Assemblage Nos. 1-4 recal- 
culated on the assumption of perfect 
mixing. Regardless of the form of the 
functional relationship between j, or 
ja and Re, comparison of Figure 17 
with Figure 6 shows that the perfect 
mixing assumption results in a very 
much greater scatter of the points 
than the assumption of no axial mix- 
ing, both for the assemblages taken 
individually and for the combined 
data. This tends to support the con- 
tention that no longitudinal mixing is 
closer to the truth than overall 
perfect mixing, in the present experi- 
ments. 

Recently, McHenry and Wilhelm 
(54) have demonstrated that their 
elaborate measurements of axial mix- 
ing of binary gas mixtures in a 
random bed of spheres can be pre- 
dicted on the assumption _ that 
“perfect mixing occurs over the 
length represented by one . 
diameter”. Based upon this model i 
can be shown (55) that the correction 
factor Y which must be applied to 
the logarithmic mean driving poten- 
tials, as in equations 1 and 2, is given 
by 


InR 


where x represents the number of 
perfect mixers in series and is there- 
fore equal to L/D,, while R is the 
ratio of inlet to outlet driving poten- 
tials. 

For the experiments, in 
which x 8, and R averages 8 for 
water partial pressure difference and 
10 for temperature difference, the 
corresponding values of Y are 0.88 
for mass transfer and 0.86 for heat 
transfer. This means that the reported 
values of j, and 
by average factors of 
respec ctively. It is 
results thus 


present 


j, should be raised 
1.14 and 1.16 
believed that the 
corrected would not 
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substantially change any of the con- 
clusions previously drawn, especially 
since the individual deviations from 
the average correction factors are 
within 4% for mass transfer and 5% 
for heat transfer. All the points in 
Figures 6 and 7 would thus be raised 
by a practically constant factor, and 
the resulting best lines would lie 
within the previously reported — 
ard deviations, with little change ir 
slope. 

The use of the logarithmic mean of 
terminal potentials is a procedure 
followed by practically all other in- 
vestigators of steady state particle- 
fluid mass or heat transfer in packed 
beds (4 6, 8, 9, 10, 18, 14, 15, 18, 23, 48, $9), 


the most notable exception being the 


differential heat transfer study of 
Denton et al. (%11), Thus the shallow 
bed results of Gamson et al. (*), 


partially summarized in Figure 9. 
would probably be subject to greater 
modification than the present data, 
due to axial mixing. The correction 
of such data reported in the literature 
by means of equation 7 is presently 
being initiated, in an attempt to 
reconcile some of the discrepancies 
between results of various investiga- 
tions, with the point data of Denton 


et al. (%11) serving as a_ possible 

datum. 

Nomenclature 

a = constant in equation (3), dimen- 
sionless. 

A = geometric surface area of spheres, 
sq. ft. 

b = constant in equation (4), dimen- 
sionless. 

c = concentration of diffusing com- 
ponent in fluid stream, Ib.-mole 
cu. ft. 

c* = solubility of diffusing component in 
fluid stream, lb.-mole/cu. ft. 

cp = molal heat capacity of diffusing 
vapor at constant pressure, B.t.u. 
(Ib.-mole) (°F.). 

C, = heat capacity of fluid stream at 
constant pressure, B.t.u./(Ib.)(°F.). 

D. = side or diameter of column, ft. 

D, = particle diameter, ft. 

D,. = moleculardiffusivity of vapor in gas, 
sq. ft./hr. 

f = friction factor = APD yg-/2LVonp, 
dimensionless. 

ge = conversion factor = 4.17 xX 108 
(Ib.) (ft.)/(hr. )? (Ib.-force) 

G = superficial mass velocity = Vop, 
Ib./(hr.) (sq. ft 

h = particle-fluid heat transfer coef- 
ficient, B.t.u./(hr.) (sq. ft.) (°F.). 

H = humidity, lb. vapor/Ib. dry gas. 

J = mass or heat transfer factor, dimen- 
sionless. 

ja = mass transfer factor 
Kees Min Sc )?/3, dimensionless. 

G 

ia’ = mass transfer factor evaluated on 
the assumption of perfect mixing 
in packed bed, dimensionless. 

jn = heat transfer factor 

h ’ , 
- (Pr)?/8, dimensionless. 
CpG 
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Min 


My, 


The Canadian Journal of Chemical Engineering, December, 1957 


= constant in equation (3) 


thermal conductivity of fluidstream, 
(B.t.u.) (ft.)/(hr.) (sq. ft.) (°F.). 
particle-fluid mass transfer coef- 
ficient, (Ib.-moles vapor )/(hr.) (sq. 
ft.) (atm). 


= length, height, or depth of packing, 


ft. 
nega- 
tive slope of log ja vs. log Re curve, 
dimensionless. 

molecular weight of 


Ib. /Ib.-mole 


molecular weight of diffusing vapor, 


Ib. /Ib.-mole. 
constant in equation (4) = 


average fluid 


stream, 


nega- 


tive slope of log jn vs. log Re curve, 


dimensionless. 


= rate of mass transfer = 


w, (H:—H)1) 
hr. 


My, lb.-moles vapor 


= partial pressure of diffusing vapor 


= mean 


= frictional pressure drop, 





in fluid stream, atm. 

partial pressure of non- 
diffusing gas in particle-fluid film, 
atm 

vapor pressure 
ponent at temperature t* 
total pressure, atm 


of dittusing com- 
atm, 


atm. 
Prandtl number Cy k, dimen- 


sionless. 
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q =rate of heat transfer = wiCp 
(ti—te) = NA + NCp (te-t*), 
B.t.u./hr. 

R =ratio of inlet to outlet transfer 
potentials: (t;—t*)/(te—t*) for 
heat transfer, (c*—c;)/(c*—c2) or 
(p*—p.)/(p*—p2) for mass trans- 
fer, dimensionless. 

Re = particle Reynolds number, DpG/p, 
dimensionless. 

Sc = Schmidt number, u/pD,, dimen- 
sionless. 

t = temperature of fluid stream, °F. 

t* = surface temperature of spheres 
during constant rate period of 
drying, °F. 

Vo = superficial velocity of fluid, ft.? 
(hr.) (sq. ft. empty column cross- 
section). 

w = weight rate of fluid flow, Ib./hr. 

w, = weight rate of non-diffusing gas 
flow, Ib./hr. 

x = “number of perfect mixers in series”’ 

= L/Dp, dimensionless. 

Y = axial mixing correction factor to 
logarithmic mean of terminal driv- 
ing potentials in mass or heat 
transfer rate equations, dimension- 
less. 

a = smallest fractional free area normal 
to major flow axis, dimensionless. 

€ = fractional void volume, external 
porosity, dimensionless. 

X = molal latent heat of vaporization 
of diffusing component at tempera- 
ture t*/B.t.u./lb.-mole. 

uM == viscosity of fluid stream, Ib/(hr.) 
(ft.). 

p = density of fluid stream, Ib./cu. ft. 

Subscripts 

1 = packed bed entrance. 

2 = packed bed exit. 

l.m. = logarithmic mean. 
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The Carbon Dioxide-Hydrogen 
Sulphide-Methane System 


Part I: Phase Behaviour at 100 F.’ 


D. B. ROBINSON? and J. A. BAILEY®* 


The phase behavior of the carbon dioxide-hydrogen sulphide- 
methane system has been studied at 100°F. and 600, 1200 and 1800 
lb./in.? abs. in a glass and stainless steel equilibrium cell of variable 
volume. The concept of the equilibrium ratio and the activity coefficient 
ratio is reviewed briefly. Analysis of the vapor and liquid phases in 
equilibrium at the three pressures permitted a calculation of equilibrium 
ratios and activity coefficient ratios as a function of pressure and 
composition. Activity coefficient ratios indicated the liquid solutions to 
be non-ideal under most of the conditions investigated. The presence of 
hydrogen sulphide suppresses the volatility of methane and increases the 
volatility of carbon dioxide whereas the hydrogen sulphide behaves 
essentially like an ideal component up to pressures of 1000 lb./in.?. 


stupy of the behavior of systems 
A containing more than one com- 
ponent over a range of temperature, 
pressure, and concentration is fre- 
quently of scientific and industrial in- 
terest. A comparison between actual 
behavior and behavior predicted for 
ideal solutions or by generalized cor- 
relations can ‘be made when experi- 
mental data exist. A knowledge of 
equilibrium phase compositions is a 
prerequisite to accurate design of in- 
dustrial equipment used to process 
multicomponent systems. 


A search of the literature failed to 
reveal any data on the ternary sys- 
tem carbon dioxide-hydrogen  sul- 
phide-methane, although the binary 
systems have been studied. It was the 
purpose of this research to investigate 
the behavior of this system at 100°F. 
and at three different pressures using 
mixtures giving a two phase system 
at the chosen conditions. Measure- 
nents made under these conditions 
permit an estimation of the nearness 
to ideality of the liquid solutions, and 


Sampling 
and 
Charging Point 


a determination of the equilibrium 
ratios for the components in the 
vapor and liquid phases. 


THEORY 
The phase behavior of single com- 
ponent, binary, and multi-component 
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systems has been reported by many 
workers. It is possible to represent 
the boundary between two phases for 
a binary system at constant tempera- 
ture as a line.on a pressure-composi- 
tion diagram. The configuration of 
this line will depend on the properties 
of the pure components at the chosen 
temperature and on the nature of the 
solutions they form. 

The addition of a third component 
to a binary system adds one degree 


Heise Gauges 







Vacuum 


Make-up Cell 


Figure 1—Schematic layout of equipment. 


of freedom. This means that the 
boundary between two phases is 
represented by a surface in three 
dimensions on a pressure composition 
digaram at a fixed temperature. It is 
not possible to represent graphically 
the variation of temperature, pressure, 
and composition simultaneously for 
a three component system. The con- 
figurations of the boundary surface 
depends on the nature of the pure 
components and the solutions they 
form. 
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Figure 2—Phase diagram at 600 lb. per sq. in. and 100°F. 


H 









x- Overall 





08 _ ~ 06 


e- Experimental Liquid 
o- Experimental Vapor 
a-Reamer, Sage and Lacey 








04 0.2 COp 


Figure 3—Phase diagram at 1200 lb. per sq. in. and 100°F. 


Data have been reported by 
Reamer and Sage ‘!) on hydrogen 
sulphide, by Sweigert, Weber and 
Allen) on carbon dioxide, and by 
Olds et al ‘*) on methane. Since the 
critical temperature for carbon di- 
oxide and methane is below 100°F., 
hydrogen sulphide is the only pure 
component existing as a liquid under 
the test conditions. 
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Bierlien and Kay ‘* report data on 
the hydrogen sulphide- -carbon  di- 
oxide system at 94 and = 105°F. 
Reamer, Sage and Lacy ‘) report 
data on the methane- hy drogen. sul- 
phide system at 100°F. The data from 
these two sources were used to con- 
struct pressure composition diagrams 
for the binary pairs at 100°F. These 
diagrams were used as a guide in 


CO, 


selecting the conditions for testing 
the behavior of the three component 
system. 


Activity Coefficient Ratios 

In ideal solutions the fugacity of 
any component in the solution in 
either the liquid or vapor phase is 


_ proportional to the mole fraction of 


the component in the phase. This is 
frequently referred to as the Lewis 
and Randall Rule ‘®. If conditions 
are such that the vapor above these 
solutions obeys the perfect gas law 
this reduces to Roult’s Law for the 
liquid phase and Dalton’s Law for the 
vapor phase. Thus 


Pi = pixi ~ (1) 


where p, is the partial pressure of 
component (i) over a liquid solution 
of mole fraction x, and vapor pres- 
sure p*,; and 


Pi = yip (2) 


where p is the total pressure of the 
system and y; is the mole fraction of 
component (i) in the vapor phase. If 
both these laws are valid 


p*ixi = pyi (3) 


These ideal relationships may 
break down for several reasons. If 
the vapor does not behave like an 
ideal gas, the partial pressure of any 
component in it will not equal the 
pressure exerted by the component 
if it alone occupied the same volume 
at the same temperature. At pressures 
significantly above the vapor pressure 
corresponding to a given temperature 
for a pure material, the effective 
vapor pressure of the liquid may be 
considerably different from the sa- 
turation pressure. Abnormal inter- 
molecular attractions or repulsions of 
different components in solution can 
also cause actual pressures to differ 
from those predicted from ideal laws. 

The activity coefficient ratio has 
been defined as a factor which allows 
for deviations from ideal behavior of 
components in solution. Thus one can 
say 


P¥i = Y'ip*ixi (4) 


where y’; is a factor which allows for 


all deviations from ideality. An al- 
ternative form is 
fyi = yif*Lxi (5) 


where pressures have been replaced 
by fugacities. In this equation f, de- 
notes the fugacity of pure com- 
ponent (i) as a vapor at the pressure 
and temperature of the system, f*;, 
denotes the fugacity of pure com- 
ponent (i) as a liquid at the tempera- 
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ture and corresponding vapor pres- 
sure of component (i) in the system, 
andy, denotes the activity coefficient 
ratio which includes the effect of 
pressure on the fugacity of the liquid 
and the effect of intermolecular 
forcé >. 

The influence of pressure on the 
fugacity of the liquid can ordinarily 
be estimated by methods which use 
the general equations for fugacity aid 
assume a constant or average specific 
volume of the liquid over the pres- 
sure range involved. In the system 
under investigation two of the com- 
ponents do not exist as liquids, so that 
the specific volume of the liquid 
imaginary. Therefore, no attempt 
was made to correct the fugacity of 
the liquid for pressure. 

Obtaining numerical values of f, 
and f*,, for use in Equation 5 involves 
the determination of a relationship 
between fugacity, pressure, tempera- 
ture, and specific volume since by 
definition 


RT dln f = vdp (6) 


at constant temperature. Upon in- 
tegration for a single component 
from its vapor pressure to any other 
pressure at a fixed temperature this 
relationship yields 


P 
fy 1 
In f*, = RT | vdp (7) 
L 
p* 


Thus, a value of the ratio f,/f*,, 
can be obtained by graphically in- 
tegrating a specific volume-pressure 
curve for any component between 
the working pressure and the vapor 
pressure at a given temperature. 
When combined with the experimen- 
tal x and y data, this ratio can be 
used to calculate activity coefficient 
ratios by rearranging equation 5 as 
follows: 


“N= >= (8) 


Equilibrium Ratios 

Equilibrium ratios are commonly 
used to express the ratio between the 
concentration of components in vapor 
and liquid phases in equilibrium at a 
fixed temperature and pressure. They 
are defined by the equation: 

yi 


Kj = (9) 
Xi 


It is readily appreciated that these 
ratios are functions of temperature, 
pressure, and composition, Equation 
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Figure 4—Phase diagram at 1800 lb. per sq. in. and 100°F. 
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Figure 5—Equilibrium ratios for methane at 100°F. 
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Figure 6—Equilibrium ratios for carbon dioxide at 100°F. 
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Figure 7—Equilibrium ratios for hydrogen 


3 shows that an ideal equilibrium 
ratio will be given by: 


Ki = (10) 


Deviations from these ideal values 
can be expected to occur for the same 
reasons that activity coefficient ratios 
are not generally unity. In terms of 
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Sage and Lacey 
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1500 


sulphide at 100°F. 


fugacity the equilibrium ratio is given 
by: 
{* 


Ki = yi 
K yi F. 


(11) 


EXPERIMENTAL METHODS 
The equipment was designed to 
permit experimental determinations 
of the variables necessary to define 


the two phase boundary surface for 
the three component system at 
100°F. A schematic layout of the 
equipment is illustrated in Figure 
It consisted of a type 316 stainless 
steel equilibrium cell A with. glass 
back and front sections capable of 
withstanding pressures of 3000 Ib./in.2 
at 200°F. and having a volume of 110 
cc. This was mounted in a constant 
temperature air bath on trunions so 
that the whole assembly could be 
rocked to speed the attainment of 
equilibrium between phases. It was 
possible to observe the phases present 
and to measure their volumes by use 
of a cathetometer reading to the 
nearest 0.1 mm. The volume of the 
sample was altered known amounts 
by the addition or removal of mer- 
cury using the high pressure cali- 
brated hand pump B. Pressures were 
measured by calibrated Heise Bour- 
don Tube gauges C and D coverin 
the ranges 0-1000 and 0-5000 Ib/in.2 
Temperatures were measured using 
copper-constantan thermocouples and 
a Leeds Northrup semi-precision po- 
tentiometer. Pressures below’ 1090 
Ib/in.2 were determined to + 1 Ib/in2 
and above 1000 lb/in? to + 3 Ib/in.’. 
The temperature was constant to 
within + 0.1°F. during the runs. 
Auxiliary equipment necessary to add 
known amounts of the pure com- 
ponents and to take samples of the 
equilibrium vapor or liquid was pro- 
vided. Connections throughout were 
made by commercial high pressure 
valves and tubing. 

The methane used in the experi- 
ments was Phillips Petroleum Com- 
pany research grade of 99.7% purity, 
the carbon dioxide obtained locally 
was 99.7% purity and the hy drogen 
sulphide ins the Matheson Com- 


pany had the followi ing analy sis: 
hydrogen sulphide 98.8%, methyl 
chloride 0.4 c, methyl mercaptan 


0.3%, nitrogen 0.2% and carbon di- 
oxide 0.3%. In order to .feduce the 
impurities admitted, the hydrogen 
sulphide was always removed from 
the liquid phase. 

A series of phase compositions was 
determined at 100°F. for three dif- 
ferent pressures, 600, 1200, and 1800 
lb/ins.2 absolute. The approximate 
overall composition of a mixture in 
the two phase region at the given 
temperature and pressure was esti- 
mated from the known two com- 
ponent boundary curves at the same 
conditions. Known amounts of each 
pure component necessary to pre- 
pare the required mixture were then 
added to the equilibrium cell. The 
amount of each component added was 
calculated from the volume, tempera- 
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ture and pressure of the sample 
using; published compressibility data 
(1,25), All lines were then purged 
with mercury and the composite 


sample was pressurized to the re- 
quirs ‘d value. The cell was rocked 
until the pressure became steady. This 
was repeated until no further de- 
clines in pressure resulted from agita- 
tion 

The volume of the liquid and 
vapor phases and the level of the 
mercury meniscus were determined 
by use of the cathetometer and a 
datum level. Following this, connec- 
tions were made to the gas sampling 
system, the top valve on the equili- 
brium cell was cracked, and mercury 
was added isobarically as the gas es- 
caped into a glass burette at atmos- 
pheric pressure until about 100 cc. 
were collected and accurately 
measured. During this operation, the 
pressure was held constant to within 
about + 3 lb/in.? of the desired value. 
Two evacuated 5 cc. glass pipettes 
were then filled with part of this 
sample and the contents analyzed for 
hydrogen sulphide by the conven- 
tional Tutwiler technique. The re- 
mainder of the sample was transferred 
to Orsat equipment where the 
methane content was determined by 
removing the carbon dioxide and 
hydrogen sulphide with potassium 
hydroxide. Carbon dioixde was then 
calculated by difference. Duplicate 
samples of vapor were always ob- 
tained. 

Without any further agitation the 
balance of the gas phase was vented 
from the equilibrium cell by slow 


TABLE 1 
EXPERIMENTAL VAPOR AND Liguip Compositions AT 100°F 


Total System Vapor Phase Liquid Phase 





Run No: aaa Naren 2 go = a, 
HS CO, CH, H.S CO. CH, H.S CO, CH, 
PREsSURE 600 Ib./in.? 
6-1 | 0.796 0.204 | 0.0 0.608 | 0.392 0.0 0.840 | 0.160 | 0.0 
6-2 | 0.792 0.143 | 0.065 0.660 0.222 0.118 0.810 | 0.170 | 0.020 
6-3 | 0.776 0.152 0.072 0.667 0.206 0.127 0.905 | 0.069 | 0.026 
6-4 0.882 0.053 0.060 | 0.684 0.115 0.201 0.937 | 0.044 0.019 
6-5 0.790 0.159 0.051 0.644 0.259 0.097 0.891 | 0.101 0.008 
6-6 | 0.804 0.196 0.0 0.630 0.370 0.0 0.815 | 0.185 | 0.0* 
6-7 0.795 0.180 0.025 | 0.630 0.317 0.053 0.870 0.122 0.008 
6-8 | 0.888 | 0.0 0.122 | 0.710 *0.0 0.290 0.967 | 0.0 0.033 
6-9 0.855 0.039 0.106 0.690 , 0.066 0.244 0.955 0.029. 0.016 
| PreEssurE 1200 lb./in.? 
12- 1 | 0.775 | 0.0 0.255 | 0.490 | 0.0 0.510 0.891 | 0.0 0.109* 
12-2 0.732 0.052 0.216 0.470 0.095 0.435 0.843 | 0.042 0.115 
12- 3 0.666 0.148 0.185 0.445 0.201 0.354 0.747 | 0.137. 0.116 
12- 4 | 0.637 | 0.216 0.147. 0.418 | 0.278 | 0.304 | 0.683 | 0.213 | 0.104 
12- § 0.546 0.306 0.147. 0.411 | 0.332 | 0.257 0.615 | 0.284 | 0.101 
12- 6 0.532 | 0.323 0.144 0.407 | 0.347 0.246 0.598 | 0.305 | 0.097 
12- 7 0.456 0.390 0.155 | 0.368 0.422 0.210 0.552 0.366 0.082 
12- 8 0.391 | 0.478 | 0.131 | 0.348 | 0.498 0.155 0.416 | 0.458 0.066 
12- 9 0.394 | 0.525 | 0.082 | 0.314 | 0.565 | 0.121 | 0.422 | 0.520 | 0.058 
12-10 |} 0.342 0.578 0.081 | 0.316; 0.579 0.105 0.372 0.563 | 0.065 
12-11 | 0.314 0.616 0.070 
Critical System Composition For 1200 Ib./in.? abs. 
PrEssURE 1800 lb./in.? 
18-1 | 0.626'0.0 0.374! 0.522 0.0 | 0.478 0.740 0.0 | 0.260 
18-2 | 0.551 0.084 0.364 Single phase 
18-3 10.590 0.059 0.351 0.530 0.047 , 0.423 0.706 | 0.036 0.258 
18-4 0.577 0.087 0.336 Single phase 
18-5 0.611 | 0.083 0.306 Single phase 
18-6 0.611 0.059 0.329 0.520 0.079 0.401 0.655 0.050 0.285 
All compositions are in mole fractions. 
*Calculated composition. 
constant pressure addition of mer- the glass burette and analyzed. It was 


possible to calculate the liquid phase 


cury. A sample of the liquid phase 
composition from the known total 


was then expanded from the cell into 


‘TABLE 2 


CALCULATED ae ILIBRIUM RATIOS AND ACTIVITY COEFFICIENT RATIOS AT 100°F 


Vapor Composition 


Cc 
H.S CO: CH, 
(mole fraction) 
0.0 0.710 0.0 0.290 
0.2 0.690 0.063 0.247 
0.4 0.680 0.128 0.192 
0.6 0.668 0.200 0.132 
0.8 0.640 0.288 0.072 
1.0 | 0.610 0.390 0.0 
0.0 0.490 0.0 0.510 
0.2 0.468 0.105 0.427 
0.4 0.439 0.225 0.336 
0.6 0.402 0.360 0.235 
0.8 0.325 0.540 0.135 
0.0 0.519 0.0 | 0.481 | 
0.2 0.550 ‘0.090 0.360 
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Liquid Composition Equilibrium Ratio Activity Coefficient Ratio 


HS CO, CH, H.S CO, CH, HS Co, CH, 
(mole fraction) 

PRESSURE 600 Ib./in.? 
0.966 0.0 |' 0.034 0.735 | 8.55 1.008 - 1.07 
0.952 0.020 0.028 0.725 3.14 8.96 0.944 1.99 1.12 
0.930 0.050 0.020 0.731 2.56 9.60 1.002 1.619 1.20 
0.896 0.090 0.014 0.745 2.22 9.43 1.021 1.404 1.08 
0.876 0.116 0.008 0.730 2.48 9.00 1.00 1.569 1.025 
0.840 0.160 0.0 0.726 2.44 0.995 1.543 

PREsSURE 1200 lb./in.? 
0.895 0.0 0.105 0.548 - 4.85 1.007 - 1.10 
0.842 0.046 0.112 0.556 2.28 3.81 1.022 2.26 0.866 
0.735 0.148 0.117 0.597 1.52 2.87 1.097 1.50 0.652 
0.585 0.315 0.100 0.688 1.14 2.38 1.283 1.13 0.541 
0.435 0.500 0.065 0.747 1.08 2.08 1.373 1.07 0.473 

PREsSURE 1800 Ib./in.? 
0.720 0.0 0.280 0.721 | Ly abana 1.484 0.555 
0.630 0.070 0.300 0.874 1.29 1.20 1.80 1.418 0.387 








©-Calculated at !00°F 
¢-Calculated at 220°F 


4—Generalized f/p plot 
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Figure 8—Fugacity ratios for ideal solutions. 
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Figure 9—Experimental and theoretical ‘K’ values at 100°F. for C 0.4. 


cell contents and the vapor phase 
volume and its analysis, so that a 
check could be obtained on the liquid 
phase analysis. 

Similar measurements were made 
at each of the three pressures. At 
1200 and 1800 Ib/in.2 the tests were 
terminated when the critical condi- 
tion was reached. 


EXPERIMENTAL RESULTS 
The experimental data obtained at 
600, 1200, and 1800 lb/in.? are pre- 
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sented in Table 1 and are indicated 
graphically in Figures 2, 3 and 4. The 
calculated liquid composition is re- 
ported in two cases where the liquid 
sample was spoiled during the 
analysis. 

At 600 Ib/in.2 the pressure is below 
the critical pressure of the hydrogen 
sulphide-carbon dioxide and hydro- 
gen sulphide-methane two phase re- 
gion at 100°F. and the two phase re- 
gion therefore extends from boundary 
to boundary. At 1200 Ib/in.? the 


pressure is above the critical pressure 
for the hydrogen sulphide-carbon 
dioxide binary but below the critical 
pressure for ‘the hydrogen sulphide- 
methane binary. Conditions similar to 
this exist at 1800 Ib°in.2 but the two 
phase region has retreated towards 
the hydrogen sulphide-methane 
boundary. 


The overall composition for run 
No. 12-11 was at or very near the 
composition of the ternary critical 
mixture+ at this pressure. -At 1205 
Ib/in.? the mixture was in the single 
phase regiong At 1203 Ib/in.? there 
was a suggestion of two phases but 
they were separated by an indistinct 
meniscus which appeared and disap- 
peared in the fluid made violently 
turbulent by the action of numerous 
tiny bubbles circulating through the 
liquid and fine droplets forming and 
falling from the vapor. A sudden re- 
duction in pressure of about one 
pound on this system resulted in such 
an eruption that the whole system 
became opaque. This separated into 
two phases separated by a very flat 
hairline meniscus after a considerable 
period of time at 1200 lIb/in.2. No 
attempt was made to sample these 
phases. The liquid phase which was 
somewhat oily at 600 lb/in.? was ex- 
ceptionally fluid at 1200 Ib/in.2. 

The pamgnennans represented by 
runs No. 18-2, 18-4 and 18-5 in the 
1800 lb/in.2 series exhibited behavior 
similar to that just described for run 
No. 12-11. They were all very close 
to the two- phase boundary surface. 
Mixtures 18-4 and 18-5 exhibited in- 
teresting color effects, becoming dark 
yellow and reddish orange instead of 
a pale straw yellow. The color reached 
its greatest intensity in run No. 18-5. 
Mixtures 18-2 and 18-4 showed 
retrograde condensation in the two 
phase region between about 1800 and 
1730 |b/in.?. 


Data of Sage, Reamer, and Lacey 
‘5) and interpolated data from Bier- 
lien and Kay “) for the binary sys- 
tems are included on Figure 2, 3 and 
4. Agreement for the hydrogen 
sulphide-methane system is good but 
repeated experiments with the hy- 
drogen sulphide-carbon dioxide sys- 
tem indicated higher concentrations 
of carbon dioxide in both the liquid 
and vapor phases. 


Equilibrium Ratios 

Equilibrium ratios were calculated 
from the experimental x-y data shown 
in Table 1. The influence of com- 
position was accounted for by using 
a parameter C defined as the ratio of 
the mole fraction carbon dioxide to’ 
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the sum of .the mole fraction carbon 
dioxide plus methane in the vapor 
phase. The use of this type of para- 
meier has been explained by Car- 
ter ‘*), Values of x and y at ‘regular 
interv -" of the parameter were read 
from Figures 2, 3 and 4. The calcu- 
lated equilibrium ratios are shown in 
Table 2 and in Figures 5, 6 and 7. 


}oundary curves at C equal to zero 
and unity representing the binary 
hy drogen sulphide-methane and hy- 
drogen sulphide- -carbon dioxide sys- 
tems respectively were determined at 
pressures intermediate to those in- 
vestigated from published data as indi- 
cated. No border curves exist at C 
equal to zero for carbon dioxide or 
for C equal to unity for methane 
since these refer to systems contain- 
ing neither component. 

Each curve is passed through the 
value K=1 at the pressure corres- 
ponding to the ternary critical pres- 
sure for the given value of C. The 
intermediate pressures were deter- 
mined from a plot of C against P 
using the experimental values at 600, 
1200, and 1800 lb/in.? and the binary 
criticals at C=0 and C=1 from 
published data. 


Activity Coefficient Ratios 

Activity coefficient ratios were 
calculated according to equation 8. 
Values of f*,,/f, were obtained for 
each component at the three experi- 
mental pressures and 100°F. by 
graphical integration. A hypothetical 
vapor pressure for carbon dioxide 
was estimated as 1210 lb/in.? and for 
methane as 6200 lb/in.* by extrapolat- 
ing a log p*, vs 1/T plot to 100°F. 
Vapor pressure data for carbon di- 
oxide were taken from Perry ‘7) for 
methane from Olds et al (8). and for 
hydrogen sulphide from Reamer and 
Sage (1), Specific volume data for 
hydrogen sulphide vapor are em- 
pirical above 394 lb/in.*, the vapor 
pressure at 100°F. The f*,f, values 
at higher pressures were eotinesed by 
extrapolating the curve using a similar 
curve calculated at 220°F. (T,= 
1.01) and data from a generalized f/p 
plot (8) as a guide. These curves are 
shown on Figure 8. 

A comparison between K_ values 
based on Dalton’s and Raoult’s Laws, 
on the ideal solution law, and on 
experiment is given in Figure 9 for 
a value of C= 0.4. 

The calculated activity coefficient 
ratios are shown in Figures 10, 11 and 
12. It is seen that increased pressure 
reduces the volatility of the more 
volatile component and increases the 
volatility of the less volatile com- 
ponent. 
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Figure 10—Activity coefficient ratios for methane at 100°F. 





Figure 11—Activity coefficient ratios for carbon dioxide at 100°F. 


Conclusions 

The hydrogen sulphide-carbon di- 
oxide-methane system forms co-exist- 
ing liquid vapor phases at 100°F. 
extending from 394 lb/in.*, the vapor 
pressure of hydrogen sulphide, to 
1910 Ib/in.?, the critical pressure of 
the binary hydrogen sulphide- 
methane system. A locus of critical 
oints exists from the hydrogen 
sulphide-carbon dioxide binary at 0.27 
mol fraction hydrogen sulphide and 
1075 Ib/in.? to the hydrogen sulphide- 
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methane binary at 0.60 mol fraction 
hydrogen sulphide and 1910 Ib/in.*. 
The activity coefficient ratios in- 
dicate the liquid solutions to be non- 
ideal under most conditions investi- 
gated. The presence of the hydrogen 
sulphide significantly influences the 
behavior of the methane and carbon 
dioxide, reducing the volatility of the 
methane and increasing the volatility 
of the carbon dioxide at all pressures 
investigated. The hydrogen sulphide 
is essentially unaffected by the pre- 
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Figure 12—Activity coefficient ratios for hydrogen sulphide. 


sence of the other components at 
pressures below about 1000 Ib/in.°. 

Work is continuing on this system 
at temperatures below the critical for 
carbon dioxide. Some _ refinements 
have been made in the experimental 
and analytical techniques. 


Nomenclature 


C = mole fraction carbon dioxide divided 
by mole fraction carbon dioxide 
plus mole fraction methane in the 
vapor phase of the ternary system, 
dimensionless. 

f = fugacity, atm. 
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mt 


- 


Xj 


H 


fugacity of a pure component at its 
vapor pressure corresponding to a 
given temperature, atm. 

fugacity of a pure component as a 
vapor at the pressure and tempera- 
ture of the system, atm. 

pressure, atm. 

partial pressure of a component i in 
a mixture, atm. 

vapor pressure of pure component i 
in a mixture at the temperature of 
the mixture, atm. 

gas constant. 

temperature, °R. 

specific volume, ft.3/mole. 

mole fraction of component at 
equilibrium in the liquid phase, 
dimensionless. 


yi = mole fraction of component i at 
equilibrium in the vapor phase, 


dimensionless. 

Y = activity coefficient ratio, dimension- 
less. 

y' = activity coefficient ratio for com- 


ponent i based on Raoult’s and 
Dalton's Laws, dimensionless. 
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Recent Developments 


in the Extrusion of Plastics’ 


G. L. BATA? 


The present state of knowledge of the operating variables affecting 
the extrusion of plastics is reviewed: An operating diagram of practical 
utility has been described and its use is recommended in connection 
with every-day problems in an extruding plant. The question of control- 
ling quality is dealt with and methods for improving quality are suggested. 


I’ took more than 2,000 years, a 
mass of empirical observations and 
trials and errors to develop the Archi- 
medean screw of the Roman’s water 
pump into the modern adiabatic 
extruder used for the processing of 
plastics. The first solid foundation for 
the quantitative picture of to-day was 
supplied in 1822 “) in the form of a 
set of basic differential flow equations. 
Many famous names established the 
solid theoretical grounds (2-4) and it 
took the effort of half-a-dozen more 
authors (5-19) until these data became 
applicable to a practical machine. 
Carley, Jepson, Mallouk, McKelvey 
and Strub “1!-17) working in a team 
collected this information a few years 
ago. They established the theoretical 
background as it is known to-day. 
This valuable work triggered an en- 
thusiastic upsurge in engineering re- 
search and development during the 
last couple of years. This paper sum- 
marizes the information presently 
available in this field. 

This paper will be confined to a 
discussion of the operation of an 
extruder and die in a simplified way. 
The description will be centered 
around polyethylene resins, as they 
are the fastest growing plastic mater- 
ial adaptable to extrusion. The find- 
ings should be true, at least qualita- 
tively, for other thermoplastic resins 
as well. 


Theoretical Background 


The essentials of an extruder are 
a cored, threaded shaft (screw) 


rotating in a fixed cylinder (barrel) 
with a feed port at one end, an exit 
port (die) at the other. The theo- 
retical data presented here apply—for 
the sake of simplicity—to a melt 
extruder, i.e. one handling preplasti- 
cized material, considered here to be 
essentially a true Newtonian liquid. 
However, for most practical purposes 
these data are applicable to an ordi- 
nary plasticizing extruder where a 
solid, granular feed is plasticized into 
a fluid. In any case, it is the melt 
section of the extruder, which pri- 
marily governs the performance of 
the machine and the resulting picture 
is thus a practically true one. 


Considering the pressure conditions 
prevailing in an extruder (Figure 
1) “8) the resin flow can be resolved 
into a forward flow component due 
to the pumping action of the screw 
and a backward flow component due 
to the pressure gradient acting against 
this pumping action, built up by the 
resistance at the die. 


This pumping action becomes ob- 
vious if the helical screw channel and 
the barrel wall are expanded into an 
equivalent straight channel, sliding at 
an angle against a plane, as on Figure 
Be 


The driving force may be resolved 
into two vector components, one 
acting parallel and the other perpen- 
dicular to the direction of the chan- 
nel. This latter is called transverse 
flow and is mainly responsible for the 


1Manuscript received August 18, 1957. 


2Associate Director of Development, Resins, Carbide Chemicals Company, Division of Union Carbide 


Canada Limited, Montreal, Que. 
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mixing action of the extruder, con- 
tributing nothing to the output. 


The whole forward flow is created 
by the vector acting alongside the 
channel, called drag flow. This drag 
flow is the most important and the 
only positive factor contributing to 
the output of the machine. It may be 
represented as ‘!5) 


tm? D?h sin @ cos @ 


Qa = — N = aN (1) 
where 
D = diameter of screw, in., 


= depth of channel, in., 
= helical angle, degrees, 


= screw speed, revolutions per sec., 


a zor 
| 


= forward flow constant, defined by 
screw dimensions as above, 


Qa = drag flow along channel, cu. in. per 
sec. a‘ 

It may be seen in Figure 3 ‘!%) that, 
the drag flow is thé sole output of the 
extruder if there is no pressure built 
up due to the resistance of a die or 
other restrictions to the flow, like a 
screen pack. 

In the case of an extruder and a 
screw of a given design its magnitude 
depends only on the rotational speed 
of the screw (9), Should the channel 
depth be changed by replacing the 
screw with a shallower one, this dis- 
charge rate will fall proportionately 
in accordance with equation 1. 
Should the screw be changed again 
to one of larger diameter, but of the 
same channel depth, the output will 
vary as the diameter. 


The backward flow created by the 
pressure built up ahead of the screw 
may also be divided into two additive 
components: the pressure flow acting 
backward in the channel and parallel 
to it and the leakage flow between 
the barrel wall and over the top of 
the screw flights, due to the clearance 
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Figure 1—Pressure profile in an extruder. 
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Figure 2—Flow pattern in extruder channel. 


existing there. These two components 
may be represented as (18) 


7Dh* Psin?d P 

0, = af (2) 

12 wL Ke 

and 
mw? D? 63 P tang Y 

Qi = ee (3) 

10eunL be 

where 

Q, = pressure flow, cu. in. per sec., 

Q; = leakage flow, cu. in. per sec., 

» = pressure drop along screw, Ibs-force 
per sq. in., 

MQ = apparent viscosity of resin, lbs-force 
sec. per sq. in 

L = effective length of screw, in., 

6 = clearance between barrel and top of 
flights, in., 

e = width of screw thread land, in 

8 = pressure flow constant defined by 
screw dimensions as above, 

Y = leakage flow constant defined by 


screw dimensions as above. 


The apparent viscosity of the 
resin at a given temperature has to 
be determined experimentally in the 


160 


extruder itself, using the equations 
previously given ‘2") at approximately 
the same speeds and output and pres- 
sure as those to which the equation 
will be applied. h has to be replaced 
by an experimentally determined H 


effective depth for an adequate 
practical description of the flow 
phenomena. 


Equations 1, 2 and 3 may be sum- 
marized for the total output (Q) of 
the extruder in a simplified equation 
as 
Q= Qa — Qp =) Qi (4) 


~ 


or 


Q =aN — BP/p — yP/p (5) 


The last member of the second 
equation may be usually neglected in 
a well designed extruder if the clear- 
ance between the top of the flights 
and the barrel wall is very small. 
Equation 5 is represented by a 
straight line if the pressure in the 
extruder head (P) is plotted against 
the output (Q), various screw speeds 


(N) giving separate lines, as in Figure 
3. This is true only if N and P are 
mutually independent variables. It is 
possible, as may be seen there, to 
impose such a great resistance to flow 
that the output of the extruder falls 
off completely. 


A characteristic output versus pres- 
sure (Q vs. P’) may be established 
for the extruder die as well ‘!9: 37); 


) = ¢(P’/M) (6) 


and if there is no screen pack used 
and the friction drop in the die head 
is negligible, P = P’ and equation 6 
may be also plotted on the same dia- 
gram of Figure 3. Theoretical form- 
ulae are available for c for various 
basic shapes such as circular, slit or 
tubular dies, but an experimentally 
determined viscosity has to be used 
in the calculations. 


Making equation 6 equal to equa- 
tion 5, the intersection of the two 
curves is the operating point of the 
extruder. The practical applicability 
of this point will be discussed later 
in this paper. 

Based on Newton’s viscosity law 
(shear stress equals the product of 
viscosity and shear rate) and consid- 
ering the vectorial implications of 
Figures 1 and 2, the total power con- 
sumption of an extruder may be inte- 
grated into the following formula (5), 


Z = (mw? D* L/h) N2u + (a/cos?@) NP + 
(7? D2e L/6 tang) N2 yu (7) 
where 


ll 


Z = power required, Ibs-force, in. per sec., 
& = apparent viscosity, lbs-force, sec. per 
sq. in. 


The same apparent viscosity may be 
used here which was previously ob- 
tained for flow equations. 

These calculations have applica- 
bility for such plastics as nylon which 
are closer to Newtonian behaviour 
than polyethylene. However they are 
applicable to the latter too until shear 

rates do not vary. Modifications of 
equation 7 have been attempted to fit 
better the behaviour of non-Newton- 
ian plastics (38) by the use of the 
power law given in equation 13. For 
practical purposes however it is ad- 
visable to use the graphical techniques 
discussed later in the section dealing 
with power requirements. 


It is possible to calculate the energy 
balance of the extruder as (2!) 


W. + Wn = Wm + Wo + Wr = (8) 


where 


W, = energy input by screw drive (con- 
sistent units), 

Wn = energy input by heaters, 

Wm = energy increase in material, 
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W.. = energy lost to cooling water of 
screw, 

W, = energy lost through radiation and 
convection. 


W.. is thé only component which 
anuot be measured directly, but it 
may be determined from Figure 4 (?!) 
as the vertical difference between the 
curves W,, and (W,, + W, — W,). 
Between the screw speeds where 
(W,, + W, W.,,) and the W, 
curves intersect the N axis the heaters 
act only as insulators to counteract 
heat losses, whereas at screw speeds 
above the W,, intersection the heaters 
switch off and the machine acts as a 
truly adiabatic extruder. In this case 
all the heat energy required is ob- 
tained from the work of the screw. 


The design of such high-speed, 
adiabatic extruders is a_ desirable 
trend, resulting in increased output 
at a high temperature and in some 
cases improved quality, not limited 
by the capacity of the barrel heaters 
and the heat sensitiv ity of the resin. 


It has been shown (16) that a well 
defined relationship can be estab- 
lished between the performance of 
geometrically similar extruders, and 
if such variables as screw speed (N) 
and -apparent melt viscosity (u) are 
kept the same, the following equa- 
tions apply: 


Q2 = x* Qi (9) 
Zz = x° Zi (10) 
P, = Pi (11) 
where 
x = D2/Di = he/hi = 62/6; = e2/e, = 
L2/Li (12) 


“ 


x” may be called scale-up factor and 
the suffix 2 refers to the larger and 
suffix | to the smaller extruder. There 
is a good agreement between ob- 
served and calculated results for melt 
extruders. 

Liquid flow may be rheologically 
described by the empirical Ostwald- 
de Waele law for non-Newtonian 
fluids (22,23), as the n-th power of the 
shear stress in proportion to the shear 
rate, if shear rates do not vary too 
widely. 


s® = K (du/dy) (13) 

Applying equation 13 to an extruder 
slit die, it may be written as 

(P’d/t)® = K (Q/ld?) (14) 


The same shear stress is expressed 
here in terms of pressure drop (P’) 
through the die, die opening (d) and 
die land length (t). The shear rate 


is given in terms of extrusion output 
(Q), slit length (1) and die opening. 
The uniformity of the extruded film 
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may be expressed as a uniformity 
index (24) 


> . St . 

U = 0; (15) 

defined as the ratio of the local ex- 
trusion rate at the far end (Q,) and 
at the feed inlet point (Q,) of the 
extrusion die. Carley ‘?4. 35) expressed 
this uniformity index as a function of 
n, 1, d, t and the radius (R) of the 
feed channel behind the die lips, thus 
presenting a workable principle for 
the rational design and operation of 
silt dies. The constants n and K are 
specific to the resin used, but are 
dependent on temperature. They may 
be easily determined experimentally 
in a suitable extrusion plastometer (26), 


It is to be emphasized that, in the 
foregoing considerations, the viscosity 







° ya ee 
AONE -p2=c(z) 
EXTRUDER OUTPUT RATE — LU.IN. [| sée. 





(u) is considered to be greatly de- 
pendent not only on the grade and 
type of the resin used, but also on 
the temperature (T). The following 
equation has been presented (39) 
describing the temperature depen- 
dence of viscosity: 


LB = Mo. e B(To-T) (16) 


where e is the base of the Naperian 
logarithms and B is the temperature 
“coefficient” of viscosity, expressed in 
reciprocal Kelvin degrees (1/K°) 


Operating an Extruder 


The two aims in operating an 
extruder for a commercial application 
are: low cost and high quality—no 
matter what is the shape or purpose 
of the product. 

Low cost consistent with high qual- 
ity, from the point of view of the 
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Figure 3—Extruder operating characteristics. 
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Figure 6—Effect of screw speed and temperature on extruder output. 


operator, can be reached only through 
increased output, good power econ- 
omy, minimum expenditure of labor, 
low maintenance cost, minimum de- 
preciation expense and the use of 
raw materials of uniform quality. 
These conditions are very closely 
interconnected, but may be broken 
down into their basic components: 
output, power, product quality and 
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design. The first three are dealt with 
here as dependent on the basic var- 
iable of the process. Due to space 
limitations, design considerations will 
not be discussed though they are con- 
sidered to be of outstanding impor- 
tance for the maximum operational 
efficiency of the machine. 

Chiefly the major process variables 
dealt with are screw speed, pressure, 


stock temperature, heating and cool- 
ing and viscosity—resulting in output 
variations, power consumption and 
variations of quality. 

Output. An experimentally de- 
termined diagram of output versus 
pressure, screw speed, resin viscosity 
(grade and temperature) is shown 
in Figure 5 “'%), This graph shows the 
major variables in a condensed and 
simple form and may be prepared 
easily for any extruder equipped with 
a pressure gauge and stock thermo- 
couples (27-25), A diagram of this 
kind is very useful in changing ex- 
truder conditions; and extrapolations 
may be safely made by the judicious 
application of the theoretical back- 
ground. 

Figure 5 is approximately equiv- 
alent to a diagram prepared for an 1} 
in, diameter extruder using a die “A” 
of c = 1.7 X 105 and a die “B” of 
Cag x Te, 

Should the extruder be running 
without a screen pack at a screw 
speed of N = 55 r.p.m. with no die 
mounted on the head of the machine, 
its operating point will be at A irre- 
spectively of the grade or temperature 
of the resin extruded. (Q = Q, + O 
+ QO). If die “A” is mounted, and a 
DYNH-type resin of melt index of 
2.0 dg/min. is used and at a stock 
temperature of 220°C, the operating 
point will shift to B and the output 
will drop from 21 to 18 lb./hr. By 
changing to a narrower die “B” at the 
same conditions, the output will drop 
further to 16.5 Ib./hr. as indicated by 
operating point F. ; 

During these losses of output the 
pressure in the extruder head in- 
creased from 0 to 500 p.s.i. due to the 
increasing die resistance and the out- 
put drop is explained by the resulting 
backward flow (Q, + Q)). 

To consider the effect of stock 
temperature, the operating point may 
be changed from B to C by lowering 
the temperature of the resin from 
220° to 180°C. This change will have 
very little effect on the output, but 
may strongly influence the quality, 
an aspect to be discussed later. The 
lowering of the stock temperature is 
equivalent in its effect to increasing 
the apparent viscosity of the resin and 
equations 5 and 6 may be called in to 
account for the changes. Similar is 
the effect of switching to a more 
viscous resin of the DYNK-type, 
having a melt index of 0.3 dg/min., as 
indicated by a change of operating 
point C to D. At first glance it may 
seem surprising that a seven times 
more viscous resin does not cut down 
the output of the extruder. The effect 
of the viscosity on output will be al- 
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mot equalized by the opposing 
effects of equations 5 and 6; but a 
closer examination of Figure 5 will 
reveal that the pressure increase will 
create an additional power require- 
ment, if the number of screw revolu- 
tions remained the same. This aspect 
will also be discussed later in further 
detail. 

lor optimum output and quality, 
the temperature of the screw is also 
of major importance beside the actual 
stock temperature itself. Figure 6 
illustrates this (31) and the graph is 
self-explanatory. The sudden drop in 
output at a high screw temperature is 
caused by the clogging and fouling of 
the helical channel of the screw by 
prematurely fused resin close to the 
feed section. This clogging eventually 
eliminates the forward flow com- 
ponent (Q, of equation 4) com- 
pletely. Such a phenomenon should 
emphasize the need for the adequate 
cooling of the feed section of the 
barrel too. 

The influence of the feed zone on 
the output of the extruder was dis- 
cussed by Barnell and Mol (2%), They 
emphasized the importance of dy- 
namic coefficients of friction between 
the particles themselves and the wall. 

Quality Considerations. The 
quality aspect of plastics extrusion is 
a rather neglected field in the tech- 
nical literature. This may be due to 
the inherent theoretical difficulties 
connected with the interpretation of 
visco-elastic flow phenomena as re- 
lated to extrudate appearance and 
properties. 

The discussion will be limited to 
quality aspects which may be con- 
trolled between the hopper and the 
die of the extruder, to avoid the 
great number of possible variable 
combinations after the extrudate left 
the die. 

One of the most frequent quality 
complaints is the roughness of the 
surface of the extrudate. Excluding 
the possibility of resin inhomogeneity 
and volatile content, Tordella (#2) 
advanced the theory of “melt frac- 
ture”.. Melt fracture is an extrudate 
roughness transverse to the extrusion 
direction, originating at the die inlet 
and occurring at the above critical 
extrusion pressures. For the sake of 
simplicity Figure 5 may be consid- 
ered again as an excellent guide to 
correct melt fracture phenomena. 
The line M,M, should represent the 
incidence of melt fracture; melt 
fracture being more and more pro- 
nounced anywhere to the right from 
this line. If extrudate roughness oc- 
curs at operating point E, this defic- 
iency may be largely eliminated by 
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Figure 7—Sectional view of the flow pattern in an extruder. 
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Figure 8—Changes in power requirement. 


slowing down the extruder to D; or 
further improvement may be ob- 
tained by changing to a softer resin 
at C. It is not absolutely necessary to 
obtain improvements at the cost of 
output, as a change from E to G will 
preserve the same productivity with 
improved quality. When none of 
these changes is beneficial or feasible, 
a tapering of the die inlet will shift 
further the line M,M, to the right of 
the chart, liberating a larger area for 
operation free of melt fracture. The 
approximate line of incidence of melt 
fracture could be conveniently de- 
termined for any extruder and die 
combination, extending Figure 5 into 
a more reliable extrusion guide. 

Haine and Land (8!) discussed var- 
ious film quality considerations as 
dependent on extruder conditions. 
Such properties as gloss, haze, clarity, 
blocking tendency, shrinkage, tensile 
strength, surging and “nerviness“ of 
the product are empirically correlated 
with extrusion variables. 
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The transverse flow in the extruder 
barrel, as previously mentioned in 
connection with equation 1, is of 
great importance from the viewpoint 
of extrudate quality. 

Combined with the leakage and 
backward flow components, it is 
responsible for the mixing and homo- 
genizing effect obtained in the ex- 
truder. This mixing effect is the 
result of a complicated flow pattern, 
represented in Figure 7 “). Back 
pressure from the die and the screen 
pack, design variables of the screw, 
and the number of screw revolutions 
per second are the controlling factors 
in the mixing action. The viscosity of 
the resin is a complex, partly inde- 
pendent variable, greatly contributing 
to the resulting extrudate homo- 
geneity. The individual effect of these 
variables may be estimated from the 
equations presented earlier, but a de- 
tailed analysis of these data should be 
the subject of a separate investigation. 
However, product quality is not de- 
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Exploded view 


pendent on this mixing effect alone, 
and the time spent by the resin in the 
extruder can not be directly corre- 
lated to quality (), 


The important problem of film 
gauge uniformity is discussed by 
Carley ‘24 25) as referred to earlier. 
By judicious application of equations 
14 and 15 and through their exten- 
sion, the following steps are recom- 
mended to improve film uniformity 
at the die, in the given order of pre- 
ference: 

a) the feed channel before the die 
lips should be of adequately large 
diameter; 

b) the land length of the die lips is 
to be varied along the length or cir- 
cumference of the die to counteract 
friction losses along the feed channel; 

c) adjustable die opening is recom- 
mended, giving a maximum attainable 
film uniformity of = 5%; 

d) die temperature may be varied 
along the slit to counteract friction 
losses (this step is not considered 
practical ). 

Carley gives quantitative data to 
guide in the use of the design and 
operating variables involved. These 
data refer to both flat and tubular slit 
dies. However, in the case of tubular 
dies, an easier solution is to com 
pletely circumvent the problem by 
the use of bottom-fed dies ‘3); these 
completely eliminate gauge uniform 
ity problems arising at the die. This 
type of die is a definite improvement 
toward an optically more uniform 
film, as gauge variation is only the 
first step to introducing further vari 


ations in haze, gloss and clarity values 


bmpirical data were published fol 
the land lengths of pipe extrusion 
dies also ‘44, in order to obtain uni 
form wall thickness. Qualitative de 
sign considerations are also discussed 


in detail elsewhere ‘4 


164 


of a slit die. 


Power Requirements. /t bas been 
mentioned earlier, that power require- 
ments may be different even at the 
same output rate, if the pressure 
created in the extruder head is 
different due to viscosity or die resis- 
tance characteristics. Considering 
Figure 8, which is a generalized case 
of Figure 5, a change of output (AQ) 
will result from a change of pressure 
(AP). The change of power require- 
ment (lb.-force X in./sec.) may be 
expressed as the product of output 
change (cu.in./sec.) and the prevail- 
ing pressure (Ib.-force/sq. in.) 


AZ =(P + AP) AQ 
liimAZ = dZ PdQ 
APO (17) 
Q 


a | pao 


| 
~ 


° 


Accordingly Z may be expressed as 
the area between the Q-axis and the 
die characteristic curve, cut off by 
horizonta! line at the operating point. 

Using a specific example taken 
from Figure 5, the change of the 
resin from the viscous DYNK (oper- 
ating point D) to the less viscous 
DYNH (operating point C) will 
result in power economy ‘at’ essen- 
tially the same output, because the 
area ODQ,O is much larger, than 
the area OCQ,O. 


Ihe actual calculation of the in- 
tegral can be avoided by either 
cutting out the respective areas from 
heavy paper and comparing their 
weights on a laboratory balance or 
by plotting the original operating 
diagram on finely divided graph paper 
and counting the number of small 
squares in both areas. However only 
differences in 


powel requirement 


should be calculated this way, as the 
areas do not represent power losses 
due to friction, etc. 


Data obtained this way are accurate 
enough for practical operating con- 
siderations. Their main value is in 
their use to extrapolate to new oper- 
ating conditions, if a change is being 
planned in the shop. Equation 7 may 
be used to counter-check results ob- 


tained. 
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The Velocity of Fall of Circulating 


and Oscillating Liquid Drops 
Through Quiescent Liquid Phases 


A. I. JOHNSON? 


This paper reports studies of the rate of fall of liquid drops passing 
through quiescent liquid phases, with an attempt to correlate the velocity 
data with the circulation and oscillation of the drops. 

Organic liquids with physical properties ranging in density from 
1,1040 to 2.9612 gms./cc., viscosity from 0.4779 to 9.630 centipoises, 
and interfacial tension from 26.08 to 45.30 dynes/cm. were used as the 
disperse phases. Water and glycerol solution ranging in density from 
0.9970 to 1.1919 gms./ce., and viscosity from 0.8937 to 27.73 centipoises 
were used as the continuous phases. In this way the more important 
physical properties affecting the fall velocity of the drops were 
investigated. 

Drops were formed from stainless steel nozzles ranging in size from 
one-half inch inside diameter to hypodermic needle size below one- 
sixty-fourth of an inch in diameter. The circulation was studied: by 
suspending aluminium particles in the disperse phase and photographing 
the drops in motion with a high speed movie camera. By this method the 
films could be projected in slow motion or studied by frame by frame 
projection. 

The voluminous data was correlated by applying a continuous 
phase viscosity correction factor to the Hu-Kintner correlation. This 
curve can be used to determine the fall velocity for a drop of any size 
whose physical properties lie within the range investigated. A break point 
in the curve serves as a criterion for oscillation. 


deals with the extension of solid 
sphere data as applied to fluid 
spheres. Recently, the velocity of fall 
of organic liquids through water has 
been studied and a correlation ob- 
tained, this 
stricted in its application. 


I* chemical engineering there are 
many operations in which two 
phases are contacted as drops of one 
phase dispersed in the other. In 
general the two fluids have different 
densities so that there is relative mo- 
tion between the disperse and con- 
tinuous phases. Liquid-liquid extrac- 
tion is one such operation for which cated that circulation and oscillation 
the velocity data and the circulation . are important factors influencing heat 
and oscillation studies of this paper and mass transfer operation. "Since 
may be of value both for basic studies there is little information on circula- 
and for the design of process equip tion and oscillation of drops, this in- 
ment. vestigation was undertaken to study 
affecting circulation and 


correlation being re- 


Recent investigations have indi- 


Literature published in the past factors 
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and L. BRAIDA®* 


oscillation and the fall of drops 
through water and liquids other than 
water as a first step for future re- 
search on heat and mass transfer 
studies. ; 


Historical Discussion 

There is considerable information 
available on the fall ged of rigid 
spheres in fluids (!:? . The friction 
factor versus Reynolds suaeilin curve 
for solid spheres has been well es- 
tablished. The motion of liquid drops 
through a continuous phase is more 
complex than that of the solid sphere, 
as evidenced by the fact that as yet 
no satisfactory correlation has been 
proposed to include existing data on 
drop motion. 

Bond “), and Bond and Newton ©), 
have shown that for small drop radii, 
the rate of descent of a drop is very 
close to the value given by Stoke's 
law for solid spheres. After a certain 
drop radius, the rate of descent 
deviates from Stoke’s law and falls 
more closely to that predicted by 
Hadamard ‘®, and Rybezynski ‘), 
who applied a correction taking into 
account viscosity and surface tension. 

Hu and Kintner (3), studied the fall 
of 10 different organic liquid drops 
through water and presented a cor- 
relation of their data with the Weber 
group, the Reynolds number, and the 
friction factor as correlating moduli. 

Licht and Narasimhamurty ‘®, con 
ducted a similar study using some 
of the liquids employed by Hu and 
Kintner. A comparison of the data 
indicates that Licht and Narasim- 
hamurty found higher fall velocities, 
the difference amounting to roughly 
106. This may be attributed to a 
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Figure 1—Apparatus for photographing circulating and oscillating drops. 


difference either in methods of 
timing employed or in the physical 
properties of the liquids used. The 
ratio of the drag coefficient of a 
liquid drop to that of a rigid sphere of 
the same volume and density and the 
same fall velocity as the drop was 
plotted against a dimensionless group 
but failed to give satisfactory results. 
Their data was also checked against 
the correlation of Hughes and 
Gilliland “, in terms of a terminal 
velocity group, a surface tension 
group and a specific gravity group, 
but poor results were obtained. In 
developing their correlation Hughes 
and Gilliland had to assume that the 
density and viscosity ratios were un- 
important. This is reasonably true 
for liquid drops falling through gases 
but not through liquids. 

Hughes and Gilliland “® have 
published a review of drop mechanics, 
in which they attempt to describe the 
effects of acceleration, distortion, cir- 
culation and oscillation to the time- 
distance-velocity relationships for the 
drop and the continuous fluid. 


Several articles on the motion of 
gas bubbles in liquids have been 
published; of these the following may 
be mentioned. Haberman and Mor- 
ton“), report work on drag and 
shape of air bubbles in water, and 
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several liquids over a very large range 
of Reynolds numbers. Garner and 
Hammerton (!2) conducted a similar 
study using bubbles of ethylene and 
carbon dioxide as well as air. Em- 
phasis was placed on the circulation 
within the bubble in accordance with 
Hadamard’s modification of Stoke’s 
law for a fluid sphere. Peebles and 
Ciarber ‘'S) have also developed a cor- 
relation for bubble velocity in various 
liquids. 

Garner ‘'*) and Garner and Skel- 
land ‘'5) studied liquid mixing as 
affected by the internal circulation 
within droplets. Their belief is that 
circulation is largely dependent on 
diffusion occurring in the drop and 
that this takes place only above a 
fairly closely defined Reynolds num- 
ber. Garner related the deformation 
of a drop to the initiation of circula- 
tion in the drop. 


Spells "© obtained photographs of 
the streamlines inside drops moving 
through another liquid, making use 
of the Schlierin pattern produced by 
imperfectly mixed water-glycerine 
mixtures, a method first used by 
Hagerty ‘!7), : 

Savic “'5) investigated circulation 
as a factor influencing the heat trans- 
fer to a drop. In this connection he 


obtained excellent photographs of 
internal circulation and_ calculated 
the time required to establish circu- 
lation when special surface effects do 
not interfere. He also considered the 
distortion of the drop suspended in a 
gas stream and found that the calcu- 
lations agree quantitatively with ex- 
periments in a vertical wind tunnel, 

Kronig and Brink ‘'%) presented a 
solution to the problem of diffusion 
in a moving spherical drop with fully 
developed circulation. The stipula- 
tion that the Reynolds number be 
one or less seriously limits the scope 
of application of their equation. 
Korchinski (29) showed that the 
Kronig-Brink theory held for Rey- 
nolds numbers greater than one. 

Korchinski studied the effect of 
circulation on heat transfer and con- 
cluded that with ‘circulating drops, 
the effective diffusivity within the 
drop is 2.25 times the molecular dif- 
fusivity. He also measured dro 
velocities and found that for Rey- 
nold’s numbers below 200, circulation 
had no effect on the friction factor 
plot of bromobenzene compared to 
the solid sphere line. Korchinski’s 
work avoided oscillating drops. 

Recently Klee and Treybal ‘?!) 
have reported rise velocities of drops 
lighter than water. They correlated 
their results using dimensional analy- 
sis for the two regions below and 
above the peak velocities indicated 
on Figures 2, 3, and 6. 


Correlation Methods 

Fluid spheres differ from solid 
spheres in that the fluid can be set 
into motion by the viscous drag in 
the fluid medium outside the drops. 
Drops are in addition subjected to 
deformation because of the dif- 
ferences in pressure acting on various 
parts of the surface. These differences 
have been the subject of many ex- 
perimental investigations and theo- 
retical analyses. 

Dimensional analysis has been used 
by many workers to obtain a cor- 
relation between drop velocity = 
the variables concerned. In this in- 
stance the following onion 
analysis was carried out giving 


= @ (Mo, Me, Pe, Ap, 7, g, D) 


D u%p, = 


o 


bo\* (odp D\" (ur g\* (Dup. 
Mc Mc Apo Ha 


Using the drop velocity data ob- 
tained in this research, various log 
plots were considered in order to 
obtain the coefficients of a, h, k, but 
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none proved feasible. By rearranging 
the above, the Np group, used by Hu 
and Kintner ‘8), was obtained: N, = = 


op scant oll 2 
By plotting Cp)yNw,Np!5 
& MeAp F 
o 
Ne 2 teh 
against vis all data for the organic 2 
E 
Oo 


dispersed phase-water systems could 















be correlated on one line similar to > 12 
Zz 
the results of Hu and Kintner. 3 
In this Hu-Kintner correlation the @ 
various terms are a8 
2 
ome ; (Po — pc) g? 3 
Drag coefficient: Cp = S 
pc u? 
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Du? e Hu & Kinter 
Weber number: Nwe ° Licht 
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| Daw % 05.06 oO? O08 09 
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In all of these numbers D is the 
equivalent diameter, that is, the 
diameter of a sphere of the same 


Figure 2—Fall velocities of carbon tetrachloride drops in water (25°C.)— 
comparison of present data with previous work. 
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volume as the drop under study. In 

this method of plotting data the 6 20 
ordinate involves the physical pro- 2 
yerties and the square of the diameter, pe. 

I P ie SS 2 
the abscissa involves the drop velo- 5 ie 
city. 

When data for organic solvent = 
drops falling through aqueous gly- 8 12 
cerol solutions was plotted on the = 
above Hu-Kintner plot, a correlation a 

: ; = 
was oO ¢ as p } e 

as not obtained as the continuous 5 “ © Carbon teleaiteride 
phase viscosity varied. A viscosity a Citereheeee 
correction for the continuous phase eo Gremebensen 
was applied with the final result that x Ethyl! bromide 
experimental data was correlated by 4 » Meee 
a plot of 

I e O-Nitrotoluene 

CpNweNp?'5 Nre + Mixture Il 
cmemens = ia against 0 
a.4 N,o-15 0 05 06 07 08 09 
bw Drop diameter cms 
Figure 3—Fall velocities of drops of organic solvents in water (25°C.) 
. a 2 
Experimental Details 
(a) Preliminary Photographic © Carbon tetrachloride 
Study of Drop Behaviour 10) @ Cilerchensens 

In order to make a comprehensive 6 Bromobenzene 
study of the behaviour of drops » Ethyl bromide 
falling through a quiescent phase, it o Miebeneen 
was sary to photograph the drop = 4 ; 

AS Necessary ph togray droy c © Gaitieeteens 
while in motion. The apparatus used = o Mkeest « Maeeed 
is shown in Figure 1. A 16 mm. Fas- = 

i ® of} v Tetrabromoethane 
tax camera, manufactured by Wol- 9 ———_——— 
lensak Optical Company, Rochester, e ne 
N.Y., was mounted on a movable 2 vo dR° 

11 f = “fe a 
piston, operated by nitrogen from a oe 
pressure cylinder which could be bay veg 
lowered at the same rate as the oO aa 
F é ‘ as . 


falling drop. A sight wire was 
mounted on the base plate in order 
that the drop could be followed. 
There are two features of this camera 
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Reynolds number 
that should be mentioned. The speed 


. Figure 4—Drag coefficients of organic solvent drops falling through water 
can be varied from zero to seven 


(25°C)—comparison with solid sphere data. 
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Kintner plot to data of this research 
—organic solvent drops falling in 
water (25°C). 


thousand pictures per second by 
varying the voltage. There is a built 
in timing light which flashes one 
hundred and twenty times a second. 
This flashing light leaves light marks 
on the film, and thus the velocity of 
fall of the drops can be determined. 
In this research the camera was 
operated at a rate of three hundred 
frames per second by two six volt 
wet cells, since the camera does not 
function very well on alternating 
current at this low voltage. 


The lighting system presented a 
problem. The falling drop could be 
photographed very well, but the cir- 
culation path indicated by the reflec- 
tion of light of the aluminium par- 
ticles could not be seen. In the final 
arrangement a light was mounted on 
an arm extended from the camera to 
provide side illumination of the drop 
as it fell. Two other photofloods 
were placed at the back of the tank. 
Dupont 930A film was found to give 
the best results. 


The dropping tank consisted of 
four pieces of plate glass cemented 
together to make a column four 
inches square. On the side of the 
tank farthest from the camera lens, 
white lines were painted every two 
centimeters from the three to the 
twenty-five centimeter mark, and 


every five centimeters from twenty- 
five to seventy-five centimeters. A 
piece of black paper was fastened 


over this so that the markings would 
show on the film. The velocity of all 
drops photographed was determined 
by counting the light blips between 
the markings on the film. 
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Figure 6—Fall velocities of organic solvent drops in aqueous 28.5% glycerol 
solution (25°C). 


Since the lights caused heating of 
the continuous phase, the tempera- 
ture was not controlled and varied 
between 25°C. and 29°C. The drops 
photographed at 25°C. were found to 
have practically the same velocity of 
fall as those obtained by stop w atch 
measurements. 


(b) Measurement of Drop Velo- 
city 

Drop velocities were measured 
a second column in which the tem- 
perature of the continuous phase 
could be controlled. A four inch 
diameter glass cylinder was placed 
inside a six inch diameter one and 
made water tight by rubber gaskets 
and brass rods connecting two brass 
plates tightly together. The space 
between the two cylinders acted as a 
constant temperature jacket, the 
water being circulated by a pump 
connected to a constant tempera- 
ture bath maintained constantly at 
25°C. + 0.2°C. Dispersed phase was 
metered from a jacketed burette to 
the nozzle ane a piece of poly- 
thene tubing, the flow controlled by 
a screw clamp. The tips of the 
nozzles were placed 3.5 cms. below 
the top of the column and immersed 
in 2.5 cms. of the continuous phase. 
Both this column and the one in 
which drops were photographed were 
large enough to avoid wall effects. 

All timing was done with a stop- 
watch manually operated and cali- 
brated against 60 cycle time signals 
to an agreement within less than one- 
tenth of a second in 10 minutes. The 
stopwatch was started as soon as the 
drop left the nozzle and stopped 


after the drop had passed a mark 
97.5 cms. down the column. As many 
as 25 velocity determinations were 
made for a single drop size each one 
varying no more than two-tenths of 
a second, and the results averaged. 


The drop volume was obtained by 
allowing 200 drops to fall and record- 
ing the volume of the aggregate. 
This was done as many times as 
needed to give reproducible results. 

When drops were formed from the 
micrometer burette, a definite volume 
was delivered and the drop then 
knocked off the needle by gently 
taping the instrument. 


(c) Systems Used 

Since density, viscosity and surface 
tension of the disperse phase in- 
fluenced the behaviour of the drops 
and similarly the density and_vis- 
cosity of the continuous _ phase, 
various disperse and continuous 
phases were investigated. The organic 
liquids chosen were essentially in- 
soluble in the continuous phases. The 
properties of distilled water were 
taken from the International Critical 
Tables, and all apparatus used in 
determining the physical properties 
was calibrated using distilled water 
at 25°C. All physical property de- 
terminations were carried out at 
25°C. £0.2°C. Interfacial tensions 
were measured by the method of 
Andreas ‘?!), 


Tables 1, 2, and 3 summarize the 
physical properties of the systems 
used. 
tetrachloride was selected 
first in order that comparisons might 


be made with previous work, It has 


Carbon 
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a viscosity near that of water and a 
density in the middle range of the 
liquids selected. Other liquids were 
chosen accordingly so that all physical 
properties could be varied at one 
time or another. Bromobenzene has a 
density close to that of carbon 
tetrachloride, but has a slightly higher 
viscosity and lower interfacial ten- 
sion. Chlorobenzene has a_ lower 
density and interfacial tension with 
nearly the same viscosity as the 
above mentioned liquids. Liquids of 
intermediate densities were prepared 
by mixing carbon tetrachloride, tetra- 
bromoethane and_ chlorobenzene; 
these are called mixtures I and II in 
the tables and on the graphs. 


(d) Method of Drop Formation 


Drops were first formed from a 
stainless steel reservoir through a 
stainless steel valve and_ thence 
through the nozzles. This was not 
suitable because aluminium particles 
accumulated at the valve and a con- 
stant flow could not be obtained. 
This method also prevented the 
simultaneous measurement of the 
drop volume. 


Glass burettes were finally used 
for measuring the volume of drops. 
Polythene tubing was used to con- 
nect the nozzles to the burette. This 
tvpe of tubing had been used by 
previous workers ‘22: 9), since it con- 
tains no soluble material which might 
alter the physical properties of the 
disperse phase. For drops of large 
volumes, a fifty milliliter burette was 
used and for smaller drops a microm- 
eter pipette was obtained. This was 
capable of delivering accurately 
volumes as small as 0.002 mls. 


The nozzle design was similar to 
that used by Hayworth and Treybal 
(23), Each nozzle was constructed of 
stainless steel carefully machined and 
bevelled at 45° to a sharp edge. These 
nozzles ranged from one-half inch 
inside diameter to one-sixteenth of 
an inch. For smaller drops hypoder- 
mic needles were used. The oblique 
end of the needle was removed and 
the tip rounded off by careful grind- 
ing on an oil stone. 


The production of constant volume 
drops presented the greatest problem 
of all in the experimental work, as 
the rate of drop formation at the 
nozzle greatly influenced the drop 
volume and velocity of fall of the 
drops specially when continuous 
phases of over thirty per cent 
glycerol in water were used. Two 
drop formation procedures were 
adopted. The first one of these was 
used for all continuous phases, while 
the second was used only for higher 


TABLE 1 


SIGNIFICANT PHYSICAL PROPERTIES FOR ORGANIC Drops FALLING IN WATER 


: Inter- 
Viscosity | Density iy facial 
iff. Tensi 
ension 
Liquid Mo Po Ap a Np 
(cos. ) (gms/cc.) | (gms/cc.) | (dynes 
cm.) 
Carbon tetrachloride. . . 0.9291 1.5857 0.5887 45.30 2.511 x 10!° 
Chlorobenzene. . . 0.7822 1.1040 0.1070 36.81 7.413 x 10% 
Bromobenzene 1.0730 1.4900 0.4930 38.12 1.786 x 10% 
Ethyl bromide 0.4777 1.4445 0.4475 31.09 | 1.067 x 10° 
Nitrobenzene 1.7839 1.1980 0.2010 26.08 1.403 x 10'° 
O-Nitrotoluene. .. . 2.0658 1.1583 0.1613 26.98 1.936 x 10% 
Tetrabromoethane. 9.3349 2.9612 1.9642 39.06 | 4.82 x 10° 
Mixture No. I*.. 2.0902 2.0511 1.0541 41.21 
Mixture No. II**.... 1.3047 1.5072 0.5102 35.63 
Carbon tetrachloride inc***. 4.4119 1.5883 0.5913 37. 84 
*A solution of tetrabromoethane, carbon tetrachloride, and chlorobenzene. 
**A solution of carbon tetrachloride and ethyl bromide. 
***\ solution of a commercial resin, Parlon, in carbon tetrachloride. 
TABLE 2 
SIGNIFICANT PHYSICAL PROPERTIES FOR ORGANIC Drops FALLING IN 
AQUEOUS 28.5% GLYCEROL 
| | | | 
| Density Density Interfacial 
Diff. Tension 
Liquid Po Ap a Np 
(gms/cc. ) (gms/ce.) ((dynes/cm.) 
Carbon tetrachloride : 1.5857 0.5172 44.26 1.286 x 10° 
Chlorobenzene. . . 1.1040 0.0355 30.93 6.39 x 109 
Bromobenzene. . . 1.4900 0.4215 34.60 7.54 x 108 
Ethyl bromide. . . 1.4445 0.3760 27.36 4.18 x 108 
Nitrobenzene.... 1.1980 0.1295 21.61 5.98 x 108 
O-Nitrotoluene. . . 1.1583 0.0898 24.39 1.24 x 108 
Tetrabromoethane. 2.9612 1.8927 38.25 2.27 x 108 
Mixture No. I... 2.0511 C.9826 38.63 
Mixture No. II.. 1.5072 0.4387 
Carbon tetrachloride inc. . 1.5883 0.5198 


Density of 28.5% Glycerol = 1.0685 gm/cc. 


Viscosity = 1.974 cps. 


TABLE 3 
SIGNIFICANT PHYSICAL PROPERTIES FOR ORGANIC Drops FALLING IN 
Aqueous 60% GLYCEROL 


Density 


Liquid Po 

| (gm/ce.) 
Carbon tetrachloride 1.5857 
Bromobenzene. . . 1.4900 
Ethyl bromide. . . 1.4445 
O-Nitrotoluene 1.1583 
Tetrabromoethane 2.9612 
Mixture No. I 2.0511 
Mixture No. II 1.5072 
Carbon tetrachloride ine. . 1.5883 


Density 600% Glycerine = 


511 gms/ce. 
Viscosity = 2 


3 cps. 


1.1 
8 


concentrations of glycerol. The first 
method consisted of forming the 
drops at a slow counting speed, one 
drop every three seconds. By this 
method the drops could be counted 
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Density Interfacial 


Diff. Tension 

Ap o Np 
(gm/ce.) ((dynes/cm.) 
0.4346 38.29 2.88 x 10° 
0.3389 J2.83 2.17 x 10 
0.2934 25.91 1.32 x 108 
0.0072 18.75 
1.8101 36.36 5.93 x 105 
0.9000 37.78 1.336 x 108 
0.3561 38.37 1.43 x 108 
0.4372 46.30 


accurately and the timing error was 
kept to a minimum. The second 
method consisted of forming the 
drops at a higher rate, thus increas- 
ing the volume slightly and the 
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Figure 7—Application of the Hu- 
Kintner plot to data of this research 


—organic solvent drops falling in 
aqueous 28.5% glycerol (25°C). 


velocity of the drop considerably. 
This rate depended on the system 
being used. Thus for concentrations 
of glycerol above 30%, there will be 
two sets of data. 

The two phases were not mutually 
saturated with one another before 
the determinations, since the mutual 
solubilities were in all cases very low. 
No difference in results were ob- 
tained by using a continuous phase 
that had been allowed to stand in the 
column for a few days. 


Observations on Drop Motion 

Because a discussion of oscillation, 
circulation and drop behaviour of all 
systems observed would be too 
lengthy, only a few general comments 
will be made. 

The frequency and manner of 
oscillation depends on the drop 
diameter, drop velocity and physical 
properties of the systems. A small 
drop with a high velocity will oscil- 
late more quickly yon a large drop 
with less velocity. A drop of tetra- 
bromoethane of 0.395 cm. diameter 
and velocity of 25.57 cms./sec. will 
oscillate quickly so that the change in 
shape can hardly be seen. In the case 
of a drop of chlorobenzene of diam- 
eter 1.086 cms. and velocity of 12.56 
ems./sec. the oscillations are very 
slow and the distortion very erratic. 

When oscillation ceases, distortion 
is present and circulation is fully 
developed. As the drop diameter and 
drop velocity decrease, the extent of 
circulation inside the drop varies 
considerably, this again depends on 
the physical properties of the two 
phases. For some drops circulation 
was observed to be fully developed 
for a short distance whereupon the 
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circulation path moved towards the 
forward stagnation point giving rise 
to what might be termed one-half 
and one-quarter circulation. For still 
smaller drop diameters, there was 
only a little motion near the forward 
stagnation point and in some cases no 
circulation present at all. 

The exact explanation for the 
decay of circulation is not yet fully 
established, as the flow patterns 
around the drop would have to be 
investigated. It might be suggested 
that there is a separation point with 
liquid spheres much like that found 
with solid spheres and that the in- 
duced circulation is greater at the 
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Figure 8— Application of the Hu- 

Kintner plot to data of this research 

—organic solvent drops falling in 
aqueous 60% glycerol (25°C). 


forward end of the drop. 

If the rate of formation at the 
nozzle was increased, the drops were 
found to have an additional velocity 
of circulation. If the drop diameter 
and velocity are sufficiently great, the 
drop will propel itself or swim down 
through the continuous phase giving 
rise to an increased velocity of fall 
and velocity of circulation even 
though the drop diameter is the same 
and circulation fully developed. This 
occurred in continuous phases having 
a concentration of glycerine greater 
than 30% 


Discussion of Graphs 
Organic Liquids Falling in Water 
Figure 2 illustrates velocity of fall 
data for carbon tetrachloride drops 
in water and shows that the data ob- 
tained lie between the results reported 
by previous workers, This difference 
is probably due to the timing methods 
employed and is small. The fall 
velocities of liquid drops increase 


linearly with diameter, reach a maxi- 
mum, and fall off asy mptotically to 
some constant value. This behaviour 
of liquid drops is different than that 
of rigid spheres. 

Figure 3 indicates similar curves 
for other liquids falling through 
water. 

The peak velocity in each case 
occurs at a drop ‘diameter where 
oscillation was observed to begin 
both visually and by photographic 
studies. The velocity plot for tetra- 
bromoethane is not shown as it is 
beyond the scope of the graph. 
Chlorobenzene is the lowest curve 
since it has the lowest density. (See 
Table 1 for physical properties.) The 
upper curves are those liquids with 
greater density differences and higher 
surface tensions. From the curves it 
appears that the effect of density on 
the fall velocity is greater than the 
effect of surface tension. V iscosity of 
the disperse phase does not appear to 
have any discernable effect on the 
fall velocity at least over the viscosity 
ranges of the dispersed phase used. 
This is evidenced by the proximity 
of the bromobenzene and ethyl bro- 
mide curves. The difference in vis- 
cosity of the two liquids is consider- 
ably greater than the difference in 
interfacial tension and density. 

The general form of the log Cp vs. 
log Nae curve is shown as Figure 4. 
At low Reynolds numbers Cp values 
fall rather closely to the curve for 
solid spheres. The slight deviation 
from the solid sphere line can be 
explained by the distortion of the 
drops. Further increase in drop size 
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brin; s about a zig zag motion and 
osciliation and deformation are pres- 
ent. Che oscillation induced by the 
peak ——r causes a sharp increase 


in the value of Cp, thus causing a 
es: rise in the log f, vs. log Re 
curve. After this abrupt rise, defor- 


mation of the drops becomes more 
and more severe and the oscillation 
progressively erratic. There is some 
crossing of the curves but all have the 
same gener ral pattern, the more dense 
liquids lying to the right while the 
lighter liquids such as chlorobenzene 
are on the left. 

The Hu-Kintner plot of Cp. 
Nye-Np®!5 against Npe/Np®> for 
this data is shown in Figure 5. A 
a break occurs at a value of 

Npe/P%!5 = 22. Again this point 
occurs at the same conditions dis- 
cussed above, that is, all drops above 
this break oscillate. 


Organic Liquids Falling in 
28.5% Glycerol 

The behaviour of drops in 28.5% 
glycerol is similar to their behaviour 
in water and velocity data is shown 
in Figure 6. The velocity increases 
linearly with drop diameter until the 
peak velocity is reached. After this 
point the velocity again decreases. 
Since the continuous phase viscosity 
has been increased the oscillations are 
not as erratic and deformation is 
quite pronounced. The curves in this 
region are considerably smoother 
than the same curves when water is 
used as the continuous phase. Chloro- 
benzene is not shown on this plot as 
the values are beyond the scope of 
the graph. As before, the upper 
curves are those liquids which have 
greater density differences, while the 
lower ones are liquids such as o-nitro- 
toluene and bromobenzene with low 
density differences. All data were 
plotted on the Cp.Ny,-P%!5 vs. 
Nae/Np®!5 graph and a mean line 
drawn through the points. This is 
shown in Figure 7. The line shows 
a definite break in the same region as 
the organic-water plot. Above this 
value all drops oscillate. By applying 
a continuous phase viscosity correc- 
tion resulting in a plot of 


(Co. Nae: Nu®!8)/(e/ bn 


the data can be made to fall on the 
water line. By this method a correla- 
tion for any system can be obtained. 


yo.t ivs, Nre N p?-15 


Organic Liquids Falling in 
60% Glycerol 

With the more viscous continuous 
phase some difficulty was encount- 
ered in obtaining reproducible results. 
With carbon tetrachloride and the 
more dense liquids fall velocity data 
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Figure 10—Modified Hu-Kintner plot for data of this research—organic solvent 
and mercury drops falling in water and in aqueous glycerol solutions. 


for small drops knocked off the mi- 
crometer pipette were not consistent 
with velocity data for drops formed 
from the other nozzles. Accordingly 
carbon tetrachloride data are not 
shown on the graphs. 


As indicated on Figure 8 data for 
other organic solvents were consis- 
tent. 

As would be expected there are no 
values at the upper end = the log 
Cp.NweNp?!5 vs. log Na. /Np0t5 
plot since there is no oscillation of 
drops. The data is considerably scat- 
tered in the upper regions of the 
graph, some point lying to the right 
of the water line. A mean line was 
drawn by eye through all points. 
Again by applying the viscosity cor- 
rection indicated the data 
coincided reasonably with the 
water line. 


above 
well 


Organic Liquids Falling in 
75% Glycerol 

In order to study 
still further 75% 
as the continuous phase. Similar re- 


the correlation 
glycerol was used 


sults were obtained and are shown in 
Figure 9. 

Korchinski’s data (25) was also used 
in checking the correlation. The only 
system that could be correlated using 
the viscosity correction was mercury 
—87.3% glycerol. Corrections to sys- 
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tems mercury — 43%, 48.8%, 97.8% 
glycerol were either above or below 
the water reference line. Since 
Korchinski did not report the densi- 
ties and viscosities of the continuous 
phases the reported results had to be 
interpolated to obtain an average 
value of densities and viscosities of 
the glycerol used. This could prob- 
ably give rise to some error. The 
more serious limitation of his data 
was that sufficient points were not 
investigated for each phase. Korchin- 
ski drew a smooth curve through all 
his data on a Hu-Kintner plot, having 
insufficient data to detect the parallel 
curves found in this research. 


Modified Hu-Kintner Correlation 
It has been _— out that the 
data in Figures 7, 8, and 9 could be 
brought down to the water line by 
dividing the ordinates by (u,/uy ois 
Using alternate points obtained in 
this research, Figure 10 was drawn 
to illustrate this modification of the 
Hu-Kintner plot. The area in which 
drops were observed to oscillate is 
also shown. The transition point be- 
tween oscillating and non oscillating 
drops is the same as the criterion Hu 
and Kintner proposed for the peak 
velocity on cures such as Figure 2. 


With a large number of points 
plotted, it is more difficult to observe 
the break in the correlation when 
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oscillation begins. However this 
region appears to be above 
Nre 
— 90, 
N p25 


Concluding Remarks 

In this paper new fall velocity data 
for organic liquid drops falling in 
water and in aqueous glycerol solu- 
tions has been presented. The 
velocity data were accompanied by 
visual and photographic studies of 
circulation and oscillation of the 
drops. 

The Hu-Kintner plot has been used 
to correlate the data. A correction 
for the viscosity of the continuous 
phase has been introduced as a divisor 
for the C, Nw.Np®!5 group proposed 
by Hu and Kintner. 

The correlating curve has been 
separated into a region of oscillating 
drops and one of non_ oscillating 
drops. There appears to be a break in 
the curve at the transition in drop 
behaviour. This corresponds to the 
point of peak velocity reported by 
Hu and Kintner. 


Nomenclature 

Cp = Drag coefficient (dimen- 
sionless ). 

D = [Diameter (cms. ). 


172 


g = Acceleration due to grav- 
ity (cms. /sec.?). 
m = Mass of a solid sphere 
(gms. ) 
o%p.2 = Physical property group 
Np = — proposed by Hu and 
gu‘Ap Kintner. 
Dup 
Nre Reynolds number 
Me 
DU 2p, 
Nwe= = Weber number. 
o 
t = Time (sec 
I = Temperature (°C.). 
u = Velocity of fall (cms./sec. ). 
Po, P = Densities of dispersed and 
continuous phases respec- 
tively (gms./cm‘), 
Ap = Density difference of sys- 
tem (gms. /cm#) 
Mo, Mt = Viscosities of dispersed and 
continuous phases (gms 
(cm) (sec. ). 
Mw = Viscosity of water 
a = [nterfacial tension (dynes 
cm), 
1, h, k = Constants. 
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Principles of Aeration as Applied 
to Waste Treatment 


Aeration systems for 


bio-oxidation of 


W. WESLEY ECKENFELDER, Jr. 


industrial wastes can be 


designed employing well-known chemical engineering principles. The 
system is designed to transfer oxygen at the same rate it is being con- 
sumed in the respiration process at a specified dissolved oxygen level. 
The oxygen requirements of the system are established from relationships 
between the BOD removed in the process and the aeration tank solids. 
To supply these requirements, aeration equipment is designed from 
previously derived sulphite oxidation studies. Factors are applied to 
these data to correct for waste characteristics which will influence the 


oxygen transfer rate. 


Economic comparisons of various aeration equipment for specific 
applications can be made employing this procedure. 


N aeration processes, mass transfer 

by diffusion occurs from gas to 
liquid (absorption) or from liquid to 
gas (stripping ) ) until a dynamic equi- 
librium is established. Mass transfer 
occurs through laminar films at the 
gas and liquid interfaces and through 
the turbulent body of fluid such that 
molecular and eddy diffusion pro- 
cesses will exist in series. Each of 
these diffusional resistances will vary 
with fluid turbulence. 


It is assumed that a dynamic equi- 
librium is established at the interface 
e.g. for each molecule which diffuses 
from gas to liquid one diffuses from 
liquid to gas). Since no concentration 
buildup will occur at the interface 
itself, the diffusional resistances may 
be considered in series and the con- 
trolling rate will be that which offers 
the greatest resistance to diffusion. It 
can be established from theoretical 
considerations that in the absorption 
or stripping of slightly soluble — 
such as oxygen and carbon dioxide i 
water, the greatest resistance occurs 
in the liquid film. 


In those systems in which the 


liquid resistance will control the rate 
of transfer (i.e. oxygen transfer in 
the activated sludge process) the 
quantity of gas transferred per unit 
time is: 


ktA (C, — Cy) (1) 


dw/dt = 


Equation 1 re-expressed in concen- 
tration units becomes: 


1 


v dw/dt =dc/dt =k, A/V (C,—CL) (2) 
where kLA/V = kya 


Due to difficulties in 
interfacial areas and surface renewal 
rates in bubble aeration the overall 
coefficient k,a is experimentally de- 
termined to define the aeration pro- 
cess. k,a is a function of bubble size 
and interfacial area, volume of liquid 
under aeration, time of contact of ai 
with liquid, tank depth, agitation 
intensity and other physical and 
characteristic of 


measuring 


chemical variables 


the system. 

In bio-oxidation processes, oxygen 
dissolved in the liquid further dif 
fuses into the microbial cell where it 
is consumed in the respiration process 


1,2), ‘To account for this oxygen 
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consumption Equation 2 must be 
modified: 


dce/dt = kta (C, — C._) — k,Sa 3 


Normal biological waste treatment 
plant operation approaches a steady 
state condition where the dissolved 
oxygen level in the aeration tanks is 
maintained substantially constant and 


> 


Equation 3 may be expresséd 
kira (C, — CL.) = k,Sa + 


In sewage and waste aeration k,a 
will be influenced by several variables 
which must be considered for process 
design and operation. For a g 
application, k,a will increase -with 
increasing air rate ewing to an in- 
crease in interfacial area and fluid 
turbulence.. In most diffusion equip- 
ment, k,a is related to gas rate by the 
exponential function 


given 


kia = CG, (5) 


Ihe gas exponent n for most com- 
mercial diffusion devices varies from 
0.6-0.9 ‘8, Some of the sparger type 
aerators have exhibited an exponent 
of 1.0-1.1 over the operating range of 
the equipment. 

The percentage of oxygen 
tion will increase with increased 


absorp 


liquid submergence, approxunately 
12° per foot of rise through the 


tank liquid \). Increasing 
ature will increase k, 


Telliper 
a due to an in 
crease in the diffusivity of oxygen in 
water and a reduction tn 
k,a may be 


VISCOSITY 
trom one 
another by the 


Cony erted 
telperature to 
relationship 


k,a at 10°C. wall be 86% and at 30°C 
1146. of the rate at 20°C 


In the 
stances and surface 


presence or ryan sub 


active agents 
the transter rate of oxygen will be 








(A) Impingement aeration equipment 
Courtesy of Walker Process 


(B) Impingement aeration equipment—tank assembly 
(Courtesy of Walker Process 


water. In the 
presence of surface active agents, an 
adsorbed 


altered from that in 


forms a barrier to 
diffusion at the air-liquid interface. 
The bubble size will also be changed, 
thus influencing the A/V ratio. Most 
substances commonly encountered in 
sewage and wastes reduce the transfer 


layer 


rate. As these substances are destroyed 
by bio-oxidation, the rate will ap 
proach that of wate 
transfel 


Some ty pic al 
industrial 


Table 1 


rates im Various 


wastes are shown in 


Measurement of K,a 


Overall 


be determined by one. of 


transfer coefhcients can 
several 
procedures. These include the sul 
phite oxidation test, non-steady state 
transfel in 


deaerated water and 


steady and non-steady state transfer 
in sludge-liquid systems 

In the presence of copper or cobalt 
salts which act as catalysts, the reac 
tion between sulphite ion and oxygen 
or alt proceeds rapidly und irrever 
sibly to completion i aqueous solu 
tion. The rate of the chemical reac 
tion is considerably more rapid than 
the absorption ol oxygen 
prepared, O.2N 


in sulphate ion and 104 


Water or waste 1 


molar in 
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Figure 1 


(C) Swing diffusion equipment—tank assembly 


Equipment Company ) 


Equipment Company ) 


cupric ion (cobalt may be substituted 
for copper). Air is applied at pre- 
determined rates and samples with- 
drawn at periodic time intervals (3-20 
minutes) depending on the air flow 
and absorption rate. The rate of 
absorption is measured by determin- 
ing the difference between the un- 
oxidized sulphite ion concentration 
before and after aeration. An aliquot 
sample, depending on the unoxidized 
sulphite ion concentration, is pipetted 
into each of two flasks containing 50 
ml. of OLN todine solution, During 
the transfer, the tip of the pipette 
should be held close to the surface 
of the iodine solution to avoid aera- 
tion. The unoxidized sulphite ion is 
then determined by a back titration 
with standard thiosulphate solution to 
a starch indicator endpoint. To mini 
mize sulphite oxidation during trans 
fer, the pipettes should be rinsed with 
distilled water and flushed with nitro 


gen from 2-3 minutes prior to samp- 


ling 
Results are expressed as p.p.m. O, 
ibsorbed per hour which is k,aC 


Von-Steady State Aeration 


Overall transfer coefhcients may be 


determined by non-steady state aera 


Courtesy of Chicago Pump Company) 


(D) Swing diffusion equipment in sewage aeration 


Courtesy of Chicago Pump Company) 


tion under standard conditions, in 
which the transfer rate is computed 
from the change in oxygen concen- 
tration of the sample under aeration. 


The dissolved oxygen present in 
the sample is stripped with a nitrogen 
purge. Air is then diffused into the 
sample at a specified rate and the 
dissolved oxygen concentration meas- 
ured at selected time intervals until 
the sample is approximately 90% 


saturated. 


The transter plot may be obtained 
over a 12 minute aeration period, ap- 
plying metered air to a liter graduate 
of sample through an Aloxite diffuser. 

‘he dissolved oxygen content may be 
ineasured either with a polarograph 
or by chemical methods. 

If the transfer character- 
istics of aeration equipment in large 


OX) gen 


tanks is to be evaluated, the use of 
nitrogen to deaerate the water is im- 
practical. In these cases sodium 
sulphite with a cobalt catalyst may 
deaerate 


be used to followed by 


aeration under non steady state cCon- 


ditions. 

Phe non-steady state absorption of 
oxygen follows the exponential func- 
tion shown in Equation 2. k,a is 
computed from the slope of a plot of 
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the saturation deficit (C, 
time of aeration. 


- C,,) vs. 


Biological Oxidation Systems— 
Steady and Non-Steady State 

In many cases it is desirable to 
evaluate k,a under operating condi- 
tions in  bio-oxidation systems. If 
steady state conditions are maintained, 
k,a can be directly computed from 
Equation 4. k,a may also be deter- 
mined employing a non-steady state 
procedure, 

The air is turned off and the dis- 
solved oxygen level permitted to ap- 
proach zero through microbial respi- 
ration. Air is then admitted at the 
desired rate and samples withdrawn 
at | minute intervals for dissolved 
oxygen measurement. Sampling is 
continued until a steady state condi- 
tion is approached. ‘The oxygen 
utilization rate and saturation value 
are determined as previously de- 
scribed. k,,a is determined by a pro- 
cedure of successive approximation, 
according to Equation 3. 


Aeration Equipment 


There are three basic types of 
aeration devices commercially avail- 
able: 

(a) Porous media type orifice diffu- 
sion units. Common types are (1) 
plates or tubes constructed of silicon 
dioxide or aluminium oxide grains held 
in a porous mass with a ceramic binder 
and (2) Saran or Nylon wrapped 
tubes or bags. These units may be 
permanently placed in the bottom of 
an aeration tank or suspended from 
flexible joints along the sidewall of a 
tank. When air is diffused through 
these units, a helical or screw motion 
is imparted to the sludge-liquid mix- 
ture. 

(b) Units employing a mechanical 
or air shear such as the impingement 
or jet aerator. 

(c) Mechanical aerators which 
entrain atmospheric oxygen into the 
sludge by surface agitation or disperse 
compressed air by a shearing and 
pumping action employi ing a rotating 
turbine or agitator. In the latter unit, 
air is discharged from a large diam- 
eter orifice in the form of a pipe or 
sparge ring beneath the agitator and 
is broken up by the shearing action 
of the high speed rotating blades of 
the agitator moving through the 
liquid. For systems of low oxygen 
utilization rate, oxygen may be sup- 
plied by air self-induced from the 
negative head produced by the rotor. 
This eliminates the 
ternal blowers or compressors. Some 


necessity for ex 


typical aeration devices are shown in 
Figure 1. 


TABLE 1 


OxYGEN ABSORPTION CAPACITY OF SOME 


Waste orSubstance cu. f 


Sewage 
after oxidation 


Kraft Mill Waste 


after oxidation 


Pulp Mill Waste 


after oxidation 


Synthetic Fibre Waste 
after oxidation 


Aeration Performance 


Standard porous diffuser units are 
designed to deliver 4-8 cu. ft./min. 
per tube. The performance of aera- 
tion equipment is generally specified 
in terms of percent absorption efh- 
ciency at a specified temperature. 
Absorption efficiency is defined as 
the Ibs of oxygen absorbed divided 
by the lbs of oxygen supplied to the 
system. The absorption efficiency 
varies with the porosity of the tube 
or plate. The absorption efficiency 
for a porosity 60 Aloxite tube in sul- 
phite solution was found to vary from 
13.5% at 4 cu. ft./min. to 12.1% at 8 
cu. fe Sanies. at 20°C. and 13.5 ft. sub- 
mergence. The Chicago Pump Co. 
has reported an absorption efficiency 
of 9.7% in sew age aer ration at an air 
flow of 4 cu. ft./min. per tube and a 


dissolved oxygen content of 0.8 
p-p-m. (20°C.). An absorption efhi- 
ciency of 8.7% was found by the 


author on activated sludge aeration 
at an average air flow of 6.4 cu. ft., 
min. The tubes are placed in the 
aeration tank 2 feet above the floor 
with a minimum and = maximum 
spacing of 6 inches and 2 feet respec- 
tively. A minimum air flow of 3 
cu. ft./min. per lineal foot of tank 
must be maintained. 

The open orifice nozzle type unit 
is designed to deliver 8-16 cu. ft./min. 
per unit at a spacing of 12-30 inches. 
Sulphite tests of this unit have shown 
a range of absorption efficiency of 
6.2% at 4 to 9.8% at 16 cu. ft./min. 
(13.5 ft. submergence). In sewage 
activated sludge, an efficiency of 
7.5% (zero dissolved oxygen) at an 
ir flow of 25 cu. ft./min. per unit 
efficiency of 8°, 
flow of 15 cu. 


was obtained. An 
was found at an air 
ft./min. per unit. Recent performance 
studies of this equipment in paper mill 
range in 
absorption efficiency of 8 to Il per 


waste oxidation showed a 


cent Over an 
24.7 cu. ft./min. 


How range of 4.5 to 


per unit at 40°C. 
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[yprcaL INDUSTRIAL WasT 
Air Flow Temp Kira Waste 
min./ 1000 cu. ft Cc 
| koa water 
aeration capacity 
10 17.0 0.82 
10 17.0 0.98 
90 36.5 0.70 
90 37.0 0.79 
20 24.0 0 70 
20 24.0 0.95 
90 20.0 2.34 
90 20.0 1.34 


The impingement aerator em 
a water stream which is air 
from the aeration tanks 
shearing device for air bubbles dis- 
charged from a large orifice mounted 
on leaders at 15-24 inch centres. The 
absorption efficiency in sulphite solu- 
tion varied from 15.2% at + to 11.5 
at 12 cu. ft./min. per unit. In Kraft 
mill waste activated sludge treatment 
the ee efficiency “was 9.4% at 
11.4 cu. ft./min. per unit. After one 
year’s omncialin the use of defoamers 
and unit clogging reduced the overall 
absorption efficiency to +.0° 


ploy S 
lifted 
acting as a 


The jet aerator pumps liquid from 
the aeration tanks through a piping 
manifold. The unit aspirates and dis- 
perses atmospheric air or air from a 
blower within the ejectors capacity. 
18.3-20°. absorption in sulphite solu- 
tion has been reported. Over a liquid 
pumping range of 34-44 U.S. gal. 
min., 0.66-1.0 Ib. 0,/hr. can be ab- 
sorbed in sulphite ‘solution. Experi- 
ments on pharmaceutical waste oxi- 
dation found a transfer efficiency of 
25% 

The turbine aerator finds principal 
application in the treatment of high 
BOD industrial wastes. While the ab- 
sorption limit of most diffused aera- 
tion units is 150 lb. 0,/dav/1000 cu. 
ft., mechanical aerators can disperse 
as much as 600 Ib. 0,/day/1000 cu. ft. 
Oxygen absorption in the turbine 
aerator will be a function of the air 
flow rate and the power input to the 
Absorption efficien 
have 


rotating agitator. 
cies as high as 30°, been re 


ported in some applications. 


Aeration Design 


(a) Oxygen Uptake Rate. [he 
oxygen uptake rate in the aeration 
system will vary as the sludge passes 
through various growth stages. he 
rate can be generally 
two phases; active respiration which 
nutrients 
growing 


considered in 


occurs in the 
and is characteristic of a 


presence or 
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Figure 2—Oxygen utilization curves through aeration tank length. 














| ° ° 

| ° 
2000 |_ 

| [- CARRE Soe 

|-Impingement Aerator -I8gpm liquor flow 

Kyo Vy / | 2- * "15 gpm 

3-Diffusair Sparger 
1000 | __ 4-Aloxite Tube - Porosity 60 

r 5- Plastic Tube 

r O Diffusair Sparger - Sewage Oxidation 

lr @ Impingement-Pulp & Paper Waste. 

b Sali a - a 
500 + Dato of Dreier Ref. (8) 

Ramee 4 4 Dich enilliten dalle ete | 

2 5 10 20 30 

G,-CU.FT./MIN. 
Figure 3—Performance curves of commercial aeration equipment. 


sludge and endogenous 
which is the autooxidation of cellular 
tissue. The respiration 
rate is 5-10% of the active respiration 
rate. 

The total oxygen requirements 
can be related to the BOD removed 
through the relationship. 


Ib. O2/day = a lb. BOD removed/day + 
blb. MLVSS) (7 


volatile 


respiré ation 


endogenous 


) 


(MLVSS is the 


solids in the aeration tanks) 


suspended 


aeration 
activated 
sludge practice, the — sludge-liquid 
mixture is rolled down the tank 
length with a spiral motion imparted 
by the air and as the BOD in the 
inlet waste undergoes oxidation and 
synthesis, the oxygen utilization rate 


In the long rectangular 
tanks used in conventional 


decreases, approaching the  endo- 
genous level at the end of the aera 


tion basins. 

Depending on the concentration of 
sludge employed and the BOD of the 
waste, the uptake rate curve through 


the aeration tanks may follow one of 
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three patterns (Figure 2). Curve 1 
would be representative of a low 
BOD waste, such as sewage, in which 
there is little or no sludge growth, 
Likewise curve 3 indicates conditions 
when the concentration of BOD is 
high and the sludge solids low, rapid 
sludge growth results in an increasing 
uptake rate until the available BOD 
is consumed. When the sludge solids 
are wir there is only a slight in- 
crease in solids through growth and 
the nis uptake rate approximates 
Curve 2. 

A variable utilization rate will 
usually not be found when square or 
circular tanks are employ ed. The 
homogenizing effect of the agitation 
and aeration will tend to equalize the 
utilization rate at the mean level. 


In order to take economic advan- 
tage of the decreasing utilization rate 
through the aeration tanks, tapered 
aeration can be employed. k,a can 
be adjusted to meet the necessary de- 
mand, by reducing the number of 
individual aeration assemblies along 


the tank length or the air rate can be 
regulated at each point in the aerator 
by appropriate valving. In deriving a 
tapered aeration design, the variation 
of k,a with oxidation must be con- 


sidered. The value of k,a relative to 
water, a, will increase with oxidation. 


(b) Oxygen Saturation. Oxygen 
saturation will be influenced by par- 
tial pressure, temperature and the 
presence of dissolved and other solids. 
At atmospheric pressure and room 
temperature in wastes, saturation will 
vary usually from 90-95% of water. 

Since air is usually released at a 
12-15 foot depth in waste oxidation 
practice, its partial pressure is in- 
Huenced by the increased pressure of 
the entering air. As the air bubbles 
rise through the tank, oxygen is ab- 
sorbed, reducing the oxygen concen- 
tration in the gas phase. The effect of 
these factors must be considered 
deriving a mean saturation value, as 
follows: 


Cate .& 
C= “ce \o04 7 2 (8) 
E oxygen saturation value 
corresponding to the average partial 
pressure of oxygen in the gas stream 
entering and leaving the aerator. tae 
is the oxygen saturation value in the 
waste at atmospheric pressure. C* i 
the oxygen saturation value in pure 
water at atmospheric pressure and 
water vapour saturation. P,, is the ab- 
solute pressure in pounds per square 
inch at the depth of air release. O, is 
the percent concentration of oxygen 
in the air leaving the tank. The effect 
of water vapor pressure is neglected 
in Equation 8. While the mean 
gradient for rising air bubbles is de- 
fined by E quation 8, the gradient for 
surface aeration will be defined by 
C, at atmospheric conditions and the 
gradient during interface formation 
by the partial pressure at the diffuser 
outlet. For practical purposes, it is 
assumed these factors produce a com- 
pensatory effect. 


C.. is the 


(ec) Oxygen Transfer. It may be 
assumed that the quantity of oxygen 
transferred per unit time by an air 
diffusion unit within the critical area 
of air dispersion is a function of 


depth, temperature, gas flow and 
waste characteristics. The transfer 
characteristics of diffusion units at 


specified depths may be related by 
the absorption number 
kia Vy 
G," 
Performance data, derived from 
sulphite oxidation at 20°C. for several 
common diffusion devices, are shown 
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in Figure 3. Data from sewage and 
pulp ind ries waste oxidation are 
also shown to illustrate the deviation 
from the sulphite curves. In Figure 
3 the solid lines represent data ob- 
tained by Dreier ‘*) in sulphite solu- 
tion at a submergence depth of 13.5 
feet and a temperature of 20°C. The 
gpm liquor flow specified is U.S. 
gal. min. 

In diffused aeration systems k, 
may be varied by (a) unit spacing 
selection; (b) diffuser permeability 
or other variable selection which de- 
fines bubble size; (c) air flow. In 
turbine aerators k,a may be varied by 


either the impellor speed o yr the air 
flow. 
Design Procedure 

1. Oxygen saturation at process 


operating temperature and _ pressure 
is computed from Equation 8. 

2. The minimum operating  dis- 
solved oxygen level should be main- 
tained between 0.5 and 1.0 p.p.m. to 
insure aerobic action. At dissolved 
oxygen levels less than approximately 
0.5 p.p.m. the respiration rate of 
the sludge is retarded and the BOD 
removal efficiency impaired. 

3. The unit oxygen uptake rate for 
the various BOD removal levels is 
estimated from an experimental plot 
or computed from Equation 7. 

4. The oxygen demand distribu- 
tion will be a function of time of 
aeration. This may be obtained by a 
laboratory or pilot plant study. The 
actual demand distribution in the 
aeration tanks will depend on the 
hydraulic characteristics of the par- 
ticular tanks as designed and will 
approximate one of the curves shown 
in Figure 2. 

5. k,a for each section of the sys- 
tem can be computed from Equa- 
tion 4. 


6. Aerator performance in sul- 
phite solution at 20°C. and 13 foot 
submergence depth can be obtained 
from Figure 3 or from similarly 
derived. The operating k,a, derived 
from Equation 5, must io corrected 
for temperature if the operating tem- 
perature deviates from 20°C. and for 
the oxygen transfer coefficient, a, of 
the specific waste. 

7. From Figure 3 (or similar plot) 
the unit air flow-volume factor to 
transfer the required oxygen is se- 
lected. 

8. After selecting a tank width and 
an air flow per unit, the unit spacing 
is computed. 

9. This should establish the design 
for the most severe operating condi- 
tions. Under less severe conditions 
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Figure 4—Aeration design chart. 


(winter operation, lower BOD load- 
ings, etc.), the reduced required air 
flow can be computed in a similar 
manner. 

The design procedure may be illus- 
trated by the following example: 
Example: 


An aeration tank is to be designed 
for the biological oxidation of an in- 
dustrial waste. The BOD removal 
and oxygen utilization rate charac- 


“teristics have been established from 


pilot plant studies. The relationship 
between BOD removal and oxygen 
utilization, Equation 7, was found to 
be: 
Ibs 0,/day =0.55 lbs BOD  re- 
moved/day + 0.08 lbs MLVSS 


If under maximum operating con- 
ditions 4,160 lbs of BOD will be re- 
moved in an aeration tank of 400,000 
U.S. gal. capacity, compute the area- 
tion equipment spacing to meet this 
requirement, assuming 2,500 p-p-m. 
MLVSS. 


(a) Oxygen Requirements 
The average oxygen requirement in 
the tank is: 


Ibs. Os/day = 0.55 (4,160) 4 
0.08 (2,500) (8.34) 0.4) 
= 2,290 + 662 


2,952 Ibs. Oo/day 
The average oxygen uptake rate is: 


2952 
r= = 36.8 ppm ‘hi 
24 X 0.4 X 8.34 


The distribution of oxygen utiliza- 
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tion is defined by curve 1, Figure 2. 
The uptake rate can be expressed as 
a ratio to the mean uptake rate in 
the tank (computed from Equation 
7). It is convenient to consider the 
distribution as average values for each 
quarter section of the tank. 


lank Ratio Uptake rate; 
Length to p.p.m. hr 

Q Mean 

0-25 ‘3 1.3. 36.8 = 47.9 
25-50 1.0 1.0. 36.8 = 36.8 
50-75 0.8 0.8. 36.8 = 29.5 
75-100 0.7 0.7. 36.8 = 25.8 


(b) Caleulation of Required « 


0.60 
- 0.90 


a for raw waste = 
a for effluent 
If the return sludge recycle is 30 
percent of the waste flow the av erage 
a at the inlet end of the tank is 0.7 X 
6.0 + 0.3 X 0.9 = 0.69 
If the increase in @ through the area- 
tion tank is linear, the mean in the 
first quarter of the tank is 0.72. 


Tank Dimensions 


The rotal tank volume is 400,000 
U.S. gal. (53,300 cu. ft.). Assuming a 
depth of 15 ft. and a width of 20 ft. 
the length is 175 ft. (44 ft./quarter). 


(ce) Aeration 


(d) Oxygen Saturation 


Operating temperature = 30°C, 
f = 0.95. Assuming transfer efficiency 
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is 8%, submergence on diffusers = 13 
ie 


iat tah. 
iit 294° 42) ~ 


0.95 x 7.63 20.3 19.6 8.4 
95 X 71.6: 94 o 2 = 8.4p.p.m. 


(e) Required k, a—first quarter 


Employing an operating dissolved 
oxygen lev el of 1.0 p-p.m. 


kia at 30°C. = ( 


(8.4 — 


kia at 20°C..= kzia (30°C. 


47 X 0.88 
7 


€ 
5 


correcting for a 
g 


kta = 5.7/0.72 = 7.92 


(f) Diffuser Spacing 

The operating data for a particular 
diffuser in water at 20°C. and 13 ft. 
submergence is as follows: 


Air flow, cu. ft./min./unit kia X V: 
610 
800 
975 
1160 
1350 


Designing for 8 cu. ft./min./unit 
kya. Ve = 1350 


V. = 1350/7.92 = 170 cu. ft. 
170 


~ = 0.565 ft. or 
15 X 20 


unit spacing = 
6.8 in. 
use average spacing of 6 inches 


The other sections of tank may be 
computed in similar fashion. If the 
computed spacing is beyond design 
limits, the unit depth and width must 
be adjusted to bring the unit spacing 
within the specification. The com- 
plete design is shown in Figure 4. 

It should be noted that in many 
cases the nature of the waste changes 
day by day and the physical condi- 
tion of the air dispersion device is 
altered by incrustation, etc. In these 
cases the performance characteristics 
of the air dispersed should be based 
on a statistical analysis of transfer 
rates under varying operating condi- 
tions. 


Conclusions 


An aeration design procedure is 
presented for activated sludge units 
treating industrial wastes. Consider- 
able work remains to be done to more 
firmly establish the aeration charac- 
teristics of various diffusion devices 
in water and to establish the nature 
and magnitude of deviations in trans- 
fer rate when these devices are em- 
ployed in industrial waste treatment. 


Nomenclature 


kia 


Cs 


lr 


= absorption coefficient, p.p.m./hr,/ 
p.p.m. concentration gradient 
oxygen'saturation in waste, p.p.m., 
corresponding to the average par- 
tial pressure of oxygen in the aera- 
tion tanks 
dissolved oxygen concentration in 
waste, p.p.m. 
fraction of BOD removed which is 
oxidized 
endogenous respiration, 
absolute temperature 
viscosity, cp. 
ratio of oxygen saturation in waste 
to oxygen saturation in pure water 
pressure at diffuser discharge, Ib./ 
sq. in. abs. 
oxygen conc entration of air leaving 
aeration tanks, ¢ 

= air flow to diffuser unit, 
ft./min. 
volume of aeration tank occupied 
by one diffusion unit, cu. ft. 
ratio of kia in waste to kza in pure 
water 
oxygen uptake rate, p.p.m./hr. 
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